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Sommaire 


L’existence d’ondes acoustiques dans l’une ou l’autre des régions du fluide que sépare l’onde 
de choc se traduit par l’apparition, dans la région en surpression, d’ondes acoustiques (réfléchies 
ou réfractées selon que les ondes incidentes sont situées dans la région en surpression ou dans la 
région en dépression) et d’ondes thermiques. Les caractéristiques de ces ondes sont calculées dans 
le cas ou l’onde incidente est plane, progressive et uniforme. Dans le cas de la réfraction l’onde 
acoustique réfractée peut, selon incidence, étre plane, progressive et uniforme, ou prendre la 
forme d’une «onde d’accompagnement» qui reste attachée au front de choc tout en glissant 
parallélement a ce front. Des constructions géométriques permettent dans tous les cas de déter- 
miner les caractéristiques cinématiques des ondes acoustiques réfléchies ou réfractées. Les relations 
dynamiques montrent que l’amplitude de l’onde réfléchie est toujours inférieure A celle de l’onde 
incidente. L’amplitude de l’onde réfractée, quelle qu’en soit la nature, peut, dans certains cas 


étre supérieure a celle de l’onde incidente. 


Zusammenfassung 


Schallwellen in dem durch eine StoBwelle gebildeten Uber- oder Unterdruckgebiet fiihren zum 
Auftreten von Schallwellen und thermischen Wellen im Uberdruckgebiet, wobei die Schallwellen 
reflektierte oder gebrochene sind, je nachdem, ob der Einfall aus dem Uber- bzw. Unterdruck- 
gebiet erfolgt. Die Eigenschaften dieser Wellen wurden fir den Fall einer ebenen fortschreitenden 
einfallenden Welle berechnet. Im Falle der Brechung kann die gebrochene Welle je nach Einfalls- 
richtung wieder eine ebene fortschreitende Welle sein oder aber die Form einer ,, Kopfwelle“ 
annehmen, die mit der StoBwellenfront fest verbunden ist und parallel zu dieser mitlauft. In allen 
Fallen 148t sich mit Hilfe geometrischer Konstruktionen die Bewegung der gebrochenen oder 
reflektierten Welle bestimmen. Aus den dynamischen Verhialtnissen ergibt sich, da die Ampli- 
tude der reflektierten Welle stets kleiner ist als die der einfallenden, wahrend die Amplitude der 
gebrochenen Welle, welcher Art sie auch sei, in gewissen Fallen gréBer werden kann als die der 
einfallenden. 


Summary 

Sound waves in one or other of the fluid regions on either side of a shock wave are made 
apparent in the region under excess pressure by acoustic waves (reflected or refracted according 
as the incident waves lie in the region of pressure excess or defect) and by thermal waves. The 
characteristics of these waves are calculated for a plane, progressive, uniform incident wave. 
In the refracted case, the refracted acoustic wave can, according to the incidence, be plane, 
progressive and uniform or take the form of an ‘‘onde d’accompagnement” which remains attached 
to the front of the shock while sliding parallel to it. Geometric constructions allow in each case 
the kinematic characteristics of the reflected or refracted acoustic waves to be determined. 
Dynamic relationships show that the amplitude of the reflected wave is always less than, but 
that of the refracted wave in certain cases greater than, the amplitude of the incident wave. 


1. DONNEES DE BASE 


1.0. Généralités 

1.0.0. Objet de travail 

La réflexion et la réfraction des ondes sonores 
par une onde de choc n’intéressent pas seulement 
Vacoustique. Ces phénoménes peuvent étre utili- 
sés en aérodynamique pour déceler ou mesurer 
Pévolution des chocs. 

Le présent travail est une tentative de théorie 
d’un cas relativement simple: onde de choc plane 


et uniforme séparant deux régions dans chacune 
desquelles le fluide est, aux phénoménes acousti- 
ques prés, en équilibre de pression et de tempéra- 
ture; onde acoustique incidente plane, progres- 
sive et uniforme. 

Ce probléme a déja fait Pobjet de quelques 
publications. Burcers [1] et BLokuinrzev [2], 
dune part, ont étudié analytiquement le cas de 
Vincidence normale dans un fluide sans viscosité. 
Ils montrent que, pour satisfaire aux conditions 
a la limite, il est nécessaire d’introduire, outre 
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l’onde acoustique réfléchie ou réfractée, une onde 
dite «d’entropie». Cette onde non progressive, 
comporte des fluctuations de température. C’est 
la forme asymptotique d’une onde thermique 


dans l’expression de laquelle on a supprimé les. 


termes d’amortissement. Les résultats ainsi ob- 
tenus en ce qui concerne cette onde ne sont donc 
valables qu’au voisinage de l’onde de choc. 

SAvER [3] d’autre part a indiqué les consé- 
quences que l’on peut tirer de la construction 
d’HvYcHENS, et montré l’existence d’incidences 
pour lesquelles aucune onde réfractée de type 
usuel (c’est-d-dire purement longitudinale) ne 
peut exister. 

L’objet du présent travail est d’examiner les 
points laissés dans l’ombre par les précédents 
auteurs, 4 savoir: comportement de l’onde ther- 
mique a grande distance, nature du phénoméne 
lorsqu’il n’y a pas d’onde réfractée de type longi- 
tudinal. Nous avons donc été conduit a traiter la 
question sous une forme plus générale, étude ana- 
lytique faite en tenant compte de la viscosité et de 
la conductibilité calorifique, incidence quelconque. 


1.0.1. Hypotheses de base et approximations 


a) Dans le calcul, nous considérons l’onde de 
choc comme une discontinuité: l’épaisseur et la 
durée d’évolution du choc sont, en effet, trés 
petites vis-a-vis des longueurs d’onde et des 
périodes acoustiques. 

b) Les équations régissant onde de choc en 
régime permanent restent valables en présence 
d’ondes acoustiques: la lenteur et la faible ampli- 
tude des phénoménes acoustiques vis-a-vis du 
phénoméne de choc justifient cette hypothése, 
tout au moins en premiére approximation. 

c) Dans l’expression des conditions 4 la limite, 
nous négligeons les petites oscillations du front 
d’onde de choc autour de sa position moyenne. 
C’est une approximation justifiée dans tous les 
problémes acoustiques: du fait que les déplace- 
ments acoustiques sont trés petits par rapport 
aux longueurs d’onde, les erreurs ainsi faites sont 
du second ordre}. 

Les conditions a la limite s’exprimeront alors 
sous forme ponctuelle, c’est-a-dire que les rela- 
tions qui les expriment ne contiennent que les 
valeurs des grandeurs en cause prises de part et 
d’autre de la limite (front de l’onde de choc) en 
un point de cette limite. 


1 Cette hypothése est implicitement contenue dans le 
calcul de Burcrers. Cet auteur calcule le mouvement de 
londe de choc par une méthode d’ailleurs inapplicable au 
cas des ondes obliques, puis, en éliminant dans les condi- 
tions aux limites les termes du second ordre, effectue une 
simplification qui correspond a notre hypothése. 


J. BRILLOUIN: ONDES DE CHOC 


ACUSTICA 
Vol. 5 (1955) 


1.0.2. Irréversibilité de l’onde de choc — 
Conséquences 


L’évolution d’une onde de choc est irréversible. 
Il est done a priori évident que les conditions a 
la limite ne pourront pas étre satisfaites si, dans 
les deux régions séparées par l’onde de choc, l’on 
mintroduit que des phénoménes acoustiques, 
c’est-a-dire adiabatiques. 

Par suite, les petits mouvements engendrés par 
la rencontre d’une onde de choc et d’une onde 
acoustique seront exprimés sous leur forme géné- 
rale telle qu’elle résulte des travaux de STOKES 
et KircHHoFF, c’est-a-dire en tenant compte de 
la viscosité et dela conductibilité thermique des 
gaz. 


1.0.3. Séparation des conditions a la limite 


Ces conditions étant ponctuelles on pourra les 
scinder en deux groupes. 

a) Conditions cinématiques, qui régissent 
la nature et la forme des ondes engendrées par 
réflexion ou réfraction. Ces conditions, lorsque 
toutes les ondes en cause sont planes, progres- 
sives et uniformes, peuvent étre mises sous la 
forme géométrique de la construction d’Huy- 
GHENS. Dans les autres cas, il est nécessaire de 
recourir au calcul pour déterminer la nature et 
la forme des ondes produites. C’est ce que l’on 
fait par exemple en optique pour traiter le pro- 
bléme de la réflexion totale et obtenir l’onde de 
CAUCHY. 

Les conditions cinématiques expriment simple- 
ment le fait que, sur la surface limite, l’évolution 
spatiale et temporelle des phénoménes qu’asso- 
cient la réflexion ou la réfraction est la méme. 
Pour les exprimer il n’est donc pas nécessaire 
d’utiliser les relations théoriques qui régissent 
Vonde de choc. 

b) Conditions dynamiques, qui n’inter- 
viennent que pour calculer l’amplitude des phé- 
noménes dont les conditions cinématiques ont 
déterminé la nature. 

Nous étudierons séparément ces deux groupes 
de conditions. 


1.1. Caractéristiques de ’onde de choe — 
Conditions a la limite 


1.1.0. Définitions, notations 


Pour désigner les deux régions du fluide sépa- 
rées par l’onde de choc nous éviterons d’utiliser 
les termes amont et aval qui, en aérodynamique 
comme en hydraulique, prétent 4 ambiguité. Si 
Vonde de choe est immobile (soufflerie, ressaut), 
cest la région aval qui est en surpression. Si 
Ponde se propage dans un fluide au repos (tube a 
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choc, onde de crue) c’est la région amont qui est 
en surpression. 
Nos notations seront donc: 


Eo = région en dépression, s’étendant de S vers 
les x positifs, 

E, = région en surpression, s’étendant de S vers 
les x négatifs, 

S = front (plan) de l’onde de choc séparant les 
deux régions, perpendiculaire a l’axe des x. 


- Nous utiliserons trois systémes d’axes: 


0 x9 y z immobiles par rapport au fluide dans la 
région Eo, 

0 x; y z immobiles par rapport au fluide dans la 
région Ej, 

_o X yz immobiles par rapport au front d’onde 

de choc S 

étant la vitesse de propagation du front 

S dans le fluide contenu en Eo, 

la vitesse de propagation de S dans le 

fluide contenu en Ej, 

U la différence entre la vitesse d’écoulement 
du fluide en FE; et sa vitesse d’écoulement 
en Eo, 


ay >0 


a,>0 


nous aurons: 


xj3=x9—Ut, X=xy—apt=x,;—a,t (1.1.1) 
et, ao — a, = U. (1.1.2) 
Nous posons enfin: 

Région E, Région E, 
pression Po Py 
masse spécifique 00 01 1.13 
température absolue To T; ee) 
vitesse du son * co Cl 


Nous prendrons comme variable indépendante 
définissant l’amplitude de V’onde de choc le 
rapport: 


&= PoP. (1.1.4) 


€ décroit de 1 4 0 4 mesure que Vamplitude du 
choc augmente. 

Nous utiliserons des fonctions sans dimension 
pour décrire le phénoméne. Nous les donnons ci- 
dessous avec leurs expressions en fonction de & 
tirées de la théorie classique de l’onde de choc. 

Posant 


pan (BY ie29]" aus 
R=® te aa) (1.1.6) 
y_U (u—1)(1—é) (1.1.7) 


eo [ut le(ut+ S]?’ 
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1/2 
, =2-| ute | 
= lwane| > 40>1 G18) 


Sow 1+yé 1/2 
A=) Fea A,<1. (1.1.9) 


Enfin, nous aurons 4 utiliser les fonctions sui- 
vantes, dans lesquelles interviennent les dérivées 


par rapport a &:R’ et Y’, de Ret Y: 


(u—1) [w+ (wt 2) é] 
2(ut &)[(u+ 1 E(w &]? ’ 


i Cea eet id 


EY'= 


(1.1.10) 


(1.1.11) 


ee uw—1?> R 
Soa estas oe 
w+ ly \(l+mé)(u+ &) 
1 - Dis pee 1/2 
ara-letrea] ihems 
aoe tl 
c=“ ery 
_ Mt (mt 2yé ( etl lie sige 
2 (u+ é) Pee wey * 
wt) , Vv 
ae ay (ev <1) (1.1.14) 
=,|("“ 2+ || eri) i 
fete jot & BH+1ljlLé(ut+é] — 


Les Tableaux I et II donnent les valeurs de 
ces diverses fonctions, calculées le plus souvent a 
la régle, précision qui suffit 4 notre propos actuel. 


1.1.1. Grandeurs ondulatoires 


La présence d’ondes de faible amplitude ayant 
une origine acoustique provoque dans les régions 
Eo et Ej de petites variations de la pression, de 
la masse spécifique et de la température autour 
de leurs valeurs moyennes. 

Pour les représenter nous posons: 


w = 0P/P= excés relatif de pression, 

s =0o0/0 = excés relatif de masse spé- 
cifique, ou condensation, 

6 =6T/T = excés relatif de température. 


(1.1.15) 


Dans le cas d’ondes acoustiques usuelles @, s 
et 9 sont trés petits et peuvent, dans le calcul, 
étre considérés comme des infiniment petits du 
premier ordre. 

Quant aux vitesses vibratoires, nous n’aurons 
a en considérer que la composante u selon ox. 


1.1.2. Conditions a@ la limite 

Les conditions 4 la limite doivent, selon nos 
hypothéses, exprimer que les relations (1.1.4) a 
(1.1.7), qui régissent l’onde de choc, sont con- 
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servées en présence d’oscillations d’origine acous- 
tique. 
Par différentiation de (1.1.4), (1.1.5) et (1.1.6) 


nous obtenons: 

@1 — Wy = —0E/E, 
dBi y Aire = 

91 — 09 = —2 5, 66 = 26 5 (1 — 9), (10,17) 


R R’ 


$1 so = R of =o (@1 — @o). (1.1.18) 


Pour obtenir, 4 partir de (1.1.7), une expres- 
sion correcte il faut prendre quelques précautions. 
On pourrait étre tenté — nous avons été — 
d’admettre, comme il est usuel en acoustique, 
que la présence des ondes acoustiques ne modifie 
pas les vitesses du son co et ¢}. 

Cette approximation reste correcte a l’intérieur 
des régions Eo et Ej, mais a leur frontiére elle ne 
peut étre admise: la relation (1.1.5) s’y oppose. 
On doit donc différentier (1.1.7) en admettant 
que co est variable. I] vient alors: 


dU = co 0Y+ Yeo (1.1.19) 


avec 0U=u —upo, différence entre les compo- 
santes uw et ug des vitesses vibratoires 


8YV Veer 2 (aon) 


D’autre part, co est donné par la relation 

classique: 
5 = 7 Poleo 
d’ou: 
2 dcq/co = 6 Po/Po — 6e0/e0 = oy — So- 

Comme nous le montrerons plus loin, dans 
notre probléme, il ne peut exister dans la région 
Eo que des ondes acoustiques incidentes. Le phé- 
noméne oscillant est done adiabatique, ce qui 
nous permet d’écrire so = @o/y, d’ou: 


dcg/co = wo/(u = 1). 


La relation (1.1.19) peut alors s’écrire: 


quwil sera commode d’utiliser sous la forme: 


pum =e (coup —H a). (1.1.21) 


I’, G, H étant les fonctions définies respective- 


ment par les relations (1.1.5), (1.1.13), (1.1.14). 


1.2. Ondes de faible amplitude dans les gaz 


1.2.0. Mise en équation 


La mise en équation de ce probléme résulte des 
travaux de SToKEs et KircHuHorr. On en trouve 
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un exposé dans louvrage de Lord RayLeicu 
(Theory of Sound, Vol. IT., § 247). Nous le résu- 


mons avec nos notations. 


u,v, Ww = composantes de la vitesse d’une particule, 
= masse spécifique du gaz au repos, 

s= 00/0 = condensation, 

Ie = pression du gaz au repos 

@® =6P/P = excés relatif de pression, 

qf = température absolue du gaz au repos, 

6=6T/T = excés relatif de température, 

y = rapport des chaleurs spécifiques, 

y = coefficient de viscosité cinématique, 

k = coefficient de conductibilité thermique égal, 


selon KircHHoFF, 4 5/2, 
vitesse du.son. 


Ss 
I 


Equation de continuité: 
ds/0t + Ou/dx + dv/éy + Ow/dz =. (1.2.1) 
Equations dynamiques: 


bay! 02s 
3 (Ot Ox 


(1.2.2) 


duc? G@ ( 3 
=? u , et cycl.eny,z 


Ep dea 
La dérivation logarithmique de léquation 
P=orT donne l’équation d’état: 


o=s+6 (1.2.3) 


Les mouvements étant trés petits on néglige 
les termes du second ordre. L’excés relatif de 
température total 6 est la somme de l’excés Oy 
résultant de la compression adiabatique et de 
Vexcés 6; provenant de la conductibilité. Or, 
on a: 


64 =(y—1)s, 


D’ou finalement l’équation thermique: 
0 / ot = (y —1) és/et + k A206. (1.2.4) 


Le calcul se conduit alors comme suit: on éli- 
mine @ en utilisant (1.2.3). On suppose ensuite 
que les diverses variables ne dépendent du temps 
que par le facteur e” (h pouvant étre complexe). 
Le facteur h remplace donc le signe de dérivation 
par rapport au temps. 

Les équations (1.2.1), (1.2.2) et (1.2.4) s’éeri- 


vent alors: 


00, /at = kA20. 


du/dx + dv/dy + dw/éz+hs=0, (1.2.1’) 
hu—vA?u=—0éM/éx, et cycl.eny,2, (1.2.2’) 
hO—kA20=(y—l)hs (1.2.4’) 

ou l’on a posé: 
M = (c2/y + hv/3) s + (c?/y) 0. (1.2.5) 


On élimine s entre (1.2.1’) et (1.2.4’), ce qui 


donne: (1.2.6) 
du/Ox + Ov/dy + Ow/éz—= —(hO—kA26)/(y—1) 
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et M s’écrit 


‘ee 1 hy ify hy 

= 1 2 
M sat (l+ 3 <r)? eee | 

Différentions les équations (1.2.2') respective- 
ment par rapport a x, y, z et ajoutons membre a 


membre. Tenons compte de (1.2.6) et éliminons 
M en utilisant (1.2.5’), il reste: 


20 hk | 4 hy 
2 
ae (ltt ga) On 
4 ee (1.2.7) 
Beata 2H?) 1g = 
ts sr) S40=0 


dont une solution particuliére est donnée par la 
forme: 


= BF Qi + Bo Qo, (1.2.8) 


ot B,, Bz sont des constantes, Q) et Q» des 
fonctions de x, y, z satisfaisant aux relations: 


A2Q;=m2Q1, A2Q2=mFQ2 (1.2.9) 
dans lesquelles m? et m3 sont déterminés par 


léquation: 


(=) ag a siege 
marry bag? a) 
heats (1.2.10 
$o(— +370) =0 


o étant le nombre sans dimension 
a = While: 


. 

Les composantes u,v,w de la vitesse vibra- 
toire dérivent alors, au facteur commun e”! prés, 
du potentiel Q; + Qo. Les constantes B; et Bz 
se tirent de (1.2.6) qui donne: 


(1.2.11) 


m- 
By =(y—1 : p=1,2). T2212 
La condensation s; correspondant a Q; est: 

si; = —(m;/h) yj e”. (1.2.13) 


L’excés relatif de pression «; s’obtient en com- 


binant (1.2.3) et (1.2.13) 


2 Baie 2) h 
m:; km; eda fone 


ae 1.2.14 
h km; —h ( ) 


oO, = — 


1.2.1. Simplifications, séparation des types d’onde, 
acoustiques et thermiques 

Dans le cas ot l’origine des ondes est acoustique, 
on aura h=iw—0d. Or, dans la relation (1.2.10) 
qui détermine les rapports m/h, h n’intervient 
que dans le nombre sans dimensions o = kh/c?. 
Lorsque ce nombre est trés petit, les racines de 
(1.2.10) prennent des valeurs approchées trés 
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simples. Dans lair, en c.g.s., v = 0,146 done 
k = 5v/2 = 0,365. Aux températures usuelles c 
est de ordre de 34 - 103 cm/s. Enfin, les pulsa- 
tions w, si on se limite supérieurement 4 la fré- 
quence, déja trés ultrasonore, 80000 Hz, sont in- 
férieures 4 0,5-106s-1. Dans ces conditions 
k w/e2 < 1,6-10-4. 

Nous avons donc le droit de considérer que le 
nombre o restera trés petit tant que le facteur 
d’amortissement 6 ne sera pas trés grand. 

Nous ferons donc les simplifications corres- 
pondantes, sauf a vérifier, lorsqu’il y aura lieu, 
que la condition o <1 est toujours respectée. 

Les racines de (1.2.10) prennent alors les 
valeurs approchées 


m2 ~(h/c)2, m2xyhi/k. (1.2.15), (1.2.16) 


La premiére nous donne une onde acoustique 
sous sa forme usuelle: Potentiel des vitesses. La 
seconde nous donne une onde de conductibilité 
thermique. 

Dans le cas d’ondes planes, la coor- 
donnée z étant éliminée par un choix convenable 
des axes, on pourra écrire: 


0; — etm,(ar+ By) (1.2.17) 


avec, pour satisfaire a (1.2.9): «2+ 62=1. 
On obtient alors le tableau de valeurs sui- 
vantes: 


Onde Onde 
acoustique thermique 


Q et? (ax+ By)/e 


etVyhik(ax+ py) 


s | (je Qer] —GmOe aay 
6 | —(y—1)s —s 

@ |ys 0 

ui Teas tVkh/y xs 

v tS TVkhl/y Bs 


1.3. Tables numériques des diverses fonctions 
introduites dans le calcul 


Tableau I 


1 0 0,714 1 1 

0,9 0,0756 0,726 1,046 0,956 
0,8 0,160 0,737 1,10 0,910 
0,7 0,26 0,756 1,17 0,86 
0,6 0,39 0,78 1,25 0,81 
0,5 0,53 0,81 1,36 0,76 
0,4 0,72 0,85 1,51 0,70 
0,3 0,96 0,92 1 73 0,63 
0,2 1,38 1,04 2,10 0,56 
0,1 2,18 1,35 2,95 0,48 
0,05 3,27 1,82 4,15 0,431 
0,025 4,75 3,18 5,86 0,405 
0 fo.) fore foe) 0,378 
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Tableau II 


1 if 

0,9 1,094 1,000 35 1,00035 
0,8 1,21 1,001 68 1,00165 
0,7 1,36 1,005 2 1,0044 
0,6 1,55 0,012 1,009 
0,5 1,80 1,025 1,018 
0,4 2,17 1,046 1,030 
0,3 2,74 1,090 1,054 
0,2 3,76 1,185 1,092 
0,1 6,18 1,45 1,166 
0,05 9,65 1,89 1,228 
0,025 11,42 2,56 1,270 

0 co co 1,323 


2. LES CONDITIONS CINEMATIQUES 
2.0. Position du probléme et formules 


2.0.0. Formes générales — Pulsation apparente Q 
et vitesse de phase w sur le front de l’onde 


de choc 


Considérons dans une des régions Ep ou Ky un 
champ acoustique incident sur S, que nous repré- 
senterons par le potentiel des vitesses 


O(r, y, 2, t). (2.0.1) 


® et x étant affectés de indice (0 ou 1) corres- 
pondant a la région en cause. 

Sur le front S ce champ produit un phénoméne 
qui, par rapport aux axes 0 X yz liés a S sera 
régi par la forme: 


Wy, z, t) = O(at, 7, 2, t) (2.0.2) 


ou a est affecté du méme indice que ®. 

Du fait que nos conditions aux limites sont 
ponctuelles, il résulte que, 4 tout champ associé 
au champ incident par réflexion ou réfraction, 
doit correspondre sur S un phénoméne dont la 
distribution spatiale et ’évolution dans le temps 
sont régis par la forme Y’. Cette condition d’ordre 
cinématique est indépendante des conditions 
dynamiques qui dérivent d’un choix particulier 
des équations qui relient les valeurs de 0, s, u, 0, 
prises de part et d’autre de S. 

Comme nous nous limitons au cas ot l’onde 
incidente est plane, progressive et uniforme, nous 
aurons: Os elalt—(zcosp + ysing)ic) 
y étant langle d’incidence. 

La fonction Y qui représente le phénoméne sur 
S prend alors la forme: 


ys — eittt—ulw) (2.0.4) 
avec 
0) 2) °) Sot ee Ay 
— = 1— A cosy = —— i=. =—tiC“ o>) 
w__2 _ l-Acosp _1—-A+(1+ 4)? 
c osing sin p a ae 
(2.0.6) 
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ot l’on a posé 
1 
idee 1 & 


c, A, w, p, tT doivent étre affectés de Vindice de la 
région qui contient l’onde. 


2.0.1. Condition pour qu’une onde acoustique plane, 
progressive et uniforme soit une onde inci- 
dente 


Considérons un mobile M entrainé par les 
ondes et cheminant sur un rayon sonore. Si, lors- 
que t croit, M se rapproche du plan S, l’onde est 
incidente. S’il s’en éloigne, l’onde est réfléchie ou 
réfractée. a. 

Par rapport aux axes liés au fluide dans la 
région considérée la vitesse avec laquelle M se 
déplace vers les x croissants est c cos@. 

Dans la région en dépression Eo, la distance de 
M aS est, a une constante prés (co cos po — ag) t. 

Elle décroit donc constamment quel que soit 
yo, puisque l’on a ag >co. Par suite, toute 
onde plane, progressive et uniforme 
existant dans Eo est une onde incidente. 

Dans la région en compression E la distance de 
M a S est, a une constante prés, (a; — c; cos¢) t. 

Or, a; < c,. Deux cas se présentent done: 


a) aj—c,cosg; <0 soit |p |< arecos Ay 
(2.0.7) 


ou |t1|<//(1—Ai)/(1+ Ai). 


L’onde est incidente. 


b) aj —c, cosy; >0 soit arc cos Ay < |gi|< 2 
(2.0.8) 


ou |r1|>/(—Ai) (A + Ap). 


L’onde est engendrée par réflexion ou réfrac- 
tion. 


2.0.2. Relations générales 


Les ondes planes, qu’elles soient acoustiques 
ou thermiques, sont régies par une forme: 


ett-mlax + By) (2.0.9) 
avec x2 + B2=] 
et m=+h/e ondes acoustiques, 
(2.0.10) 
m=+)/yh/k ondes thermiques. 


Sur le plan S qui se déplace, dans la région 
considérée, a la vitesse a, l’'exposant de (2.0.9) 
doit prendre la forme i (t—y/w); d’ow les 


relations: 
iQ=h—xma, iQ/w=Bm (2.0.11) 
c’est-a-dire 
a=(h—iQS“)/am, B=iQ/wm (2.0.12) 


——— 
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ax? + 62 —1 donne 
(a m/Q2)? + (a/w)2 — (h/Q —i)? = 0. (2.0.13) 
En remplagant dans cette équation m par sa 
valeur tirée de l’une ou I’autre relation (2.0.10) 
on obtient une équation qui déterminera h en 
fonction de 2 et w. On tire ensuite m, x et 8 de 


(2.0.10) et (2.0.12). 


2.1. L’onde thermique réfléchie ou réfractée 


Portant dans (2.0.13) la valeur m?=yh/k il 
vient: 


(2.1.1) 

(h/Q)2 —2 (i+ ya2/2 k Q) h/Q—(1—a2/w?) = 0 

aon: (2.1.2) 
hiptoiys? / 2kQ\? ey 
romana 2 80 ji 14 (oS) +2 ya2 |" 


Pour nous convenir, la racine doit avoir une 
partie réelle négative: onde éyanouissante lorsque 
t croit indéfiniment. 

Or, nous l’avons vu, le nombre sans dimensions 
k Q/a? est trés petit. D’autre part, (2.0.6) donne 
Q/w = (w/e)sing (w,¢, gy affectés de Vindice de 
la région qui contient l’onde incidente). 

Par suite, k Q/aw = (kw/ac) sing est du méme 
ordre de grandeur que k Q/a?, c’est-a-dire trés 
petit. 

Prenant alors la valeur approchée du radical 
de la formule (2.1.2) on voit que la racine qui 
nous convient est celle ot ce radical est affecté 
du signe négatif. 

Pour obtenir une valeur approchée de h/Q, il 
faut développer le radical jusqu’aux termes du 
second ordre, les termes du premier ordre dis- 
paraissant, on obtient alors: 


kQ2 az+ w2 


On voit aisément que 
kh 1 222 (a2 + w?) 
ce yee Roa 


reste toujours trés petit. L’approximation faite 
pour séparer l’onde thermique de l’onde acous- 
tique est donc valable. 

D’autre part (2.1.3) montre que h est réel et trés 
petit. Par suite, l’onde thermique n’est pas sen- 
siblement oscillante et son amortissement est lent. 

De la valeur de h on tire: 


m= iQVa? + w2/aw, (2.1.4) 
w 
— teat Va a = cos Pr, 
(2.1.5) 
B <= sing, 
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La distribution spatiale des excés de tempéra- 
ture et des condensations est donc sinusoidale. 
L’onde est oblique. Elle n’est pas progressive et 
s’éteint sur place. 

Il en résulte qu’elle ne peut étre engendrée que 
derriére le plan S (dont elle ne peut précéder le 
mouvement). Elle ne peut donc exister que dans 
la région en compression FE}. 

L’absence de toute onde thermique 
dans la région Ep justifie done la mé- 
thode par laquelle nous avons établi la 
condition aux limites (1.1.20). 

Le potentiel Q =e” %***"”) s’écrit alors: 


(2.1.6) 


Les grandeurs liées a l’onde sont alors (voir 


1.2.13): 


Q aes Ao (z/a — yw) 


= — 7 Qett, 6; = —s;, wo; = 0, 


et les composantes de la vitesse (2.1.7) 


Ces valeurs nous seront utiles pour la simpli- 
fication des conditions dynamiques. 


2.2. L’onde acoustique réfléchie 


Dans ce cas |’onde incidente est située, comme 
Vonde réfléchie, dans la région en surpression Ej. 
Les relations (2.0.5) et (2.0.6), qui lient les carac- 
téristiques cinématiques w et p (ou t) de l’onde 
acoustique a la pulsation © et la vitesse de 
phase w observables sur le front d’onde de choc 
sont valables 4 la fois pour l’onde incidente 
(m1, 1, T1) et Vonde réfléchie (w}, 7}, T})- 

Eliminant 2 et w on obtient: 

T71 T= (1—Aj)/(1+ Aj), 
oO} sing, _ (1— A,)?+ (1 + A?) Te 


@, sing; (1—Ay)(1+ A1)(1+ 72)" 


(2.2.1) 


(2.2.2) 


(2.2.1) montre que, a l’onde incidente (1; < 
V(1—A1)/(1+A1) correspond bien une onde 
réfléchie [1 > /(1— A1)/(1+ A)) J. 

(2.2.2) donne le rapport des pulsations (pergues 
par un observateur entrainé par le fluide) de 
l’onde incidente et de londe réfléchie. 

Lorsque ¢ croit de 0 (incidence normale) a 
arccos A, (angle d’incidence limite) 1; croit de 0 a 


V(1—Ay)/(1+ Aj), @}/e1 eroit de (1—A})/(1+A1) 
al. 
On a done toujours: 


ep a 
ii 0 Coa 
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Du fait que toute onde incidente plane, pro- 
gressive et uniforme engendre une onde réfléchie 
de méme type, la construction d’HUYGHENS est 
valable. 


2.3. L’onde acoustique réfractée 


2.3.0. Généralités — Séparation des deux types 
d’ondes réfractées 


Pour qu'il y ait réfraction il faut que l’onde 
incidente soit située dans la région en dépression 
Eo. D’autre part toute onde plane progressive 
et uniforme située dans Eo est une onde incidente 
et donnera lieu a la production d’ondes réfractées. 

La vitesse de phase w du phénoméne engendré 
sur S par l’onde incidente est donnée par (2.0.6): 


Vi) (2 a] c 
2 To 


(2.3.1) 


ou Ag>l. 


x 


Lorsque l’angle d’incidence yo varie de —z a 
+2, To varie de —oo A + co. w croit constamment 
et sa variation est donnée par le tableau suivant: 


—arecos1/Ao 0 arccos1/Ao 


yo | —7 


fore) je 0 
¥9 Ag+ 1 


| 
Ao+1 


cae RAS 


w | —oco 0 jigs) = 6S 0 


La méme vitesse de phase w peut donc étre 
produite par deux ondes dont les incidences sont 
comprises, l’une entre —z et 0, l’autre entre 0 et 
+a, Pune étant en valeur absolue supérieure 4 
arccos1/Ao, autre inférieure a cette valeur. Les 
valeurs correspondantes de T; 7 et 7, relatives 
a ces deux ondes, sont liées par la relation: 


TO ™ = —(Aop—1)/Ao+ ib (2.3.2) 


Enfin, (2.0.6) nous montre que les pulsations 
des deux ondes incidentes qui produisent sur S 
le méme phénoméne (égalité de w et de Q) sont 
dans le rapport: 


0/@o = |singo/sin ¢6). 


Ainsi, le phénoméne acoustique observé sur S, 
et a fortiori a l’intérieur de la région E, peut étre 
engendré indifféremment par l’une ou l’autre, ou 
toute combinaison linéaire, des deux ondes inci- 
dentes ci-dessus. Aucune observation faite en E; 
ne peut lever le doute. 

La nature de l’onde réfractée dépend de la 
valeur absolue de la vitesse de phase w. Si | w| > 


cj — a} Vonde réfractée est plane, progressive 
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et uniforme. Si |w| </c?—a?, on ne peut lui 

faire correspondre dans E; aucune onde plane, 

progressive et uniforme, L’onde réfractée prend 

alors une structure exponentielle et nous la 
nommerons onde d’accompagnement. 

Les angles d’incidence go correspondant a la 

. 


valeur de w pour laquelle on passe de lun a 
l’autre type d’onde réfractée sont donnés par: 


Vee + Vee —a? + a? — ce 
Re : (2.3.3) 


To = 


Lorsque |to| est extérieur a l’intervalle défini 
par les limites ci-dessus, l’onde réfractée est plane, 
progressive et uniforme. Dans le cas contraire on 
obtient onde d’accompagnement. 

Exprimées en fonction des relations théoriques 
qui régissent l’onde de choc, les deux valeurs 
limites de To s’écrivent: 


(2.3.4) 
ats ee tVI+ HE + V4 DOtS 
ahs Vuté+V(utbé 


Les valeurs des angles limites correspon- 
dants sont données dans la Fig. 1. 


\ 


Zone de Vonde 


05 i 
d’ accompagnement 


0° 45° 90° 135° 
——> Angle d’ incidence op, 


Fig. 1. 


2.3.1. Onde réfractée plane, progressive et uni- 
forme: w? > ci —a? 
L’équation (2.0.6) appliquée aux deux régions 
Eo et E, qui contiennent respectivement l’onde 
incidente et l’onde réfractée donne la relation qui 


lie to et T: 
(2.3.5) 


co— ao + (co + ao) t 


e1—ay + (e1 +. a1) ty 


T1 TO 
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to étant imposé par Vincidence qo, tT) est déter- n|Q = —A /(1— A?) =—Aj@4/Q, (2.3.14) 
miné par une équation du second degré dont : 18 
seule la racine supérieure A |/(cy—aj) /(c1 + a1) a shane (2.3.15) 
correspond a l’onde réfractée. L’autre racine ‘ 
donne dans Ej l’onde incidente qui pourrait par ny’ /|Q = —c\/v, (2.3.16) 
réflexion engendrer la premiére. me 
On a donc: ae A) 
C9 — a9 + (¢9 + a0) Te + V[co—ao + (co + ao) t)]2—4 (c? — a?) ot 
0s Ma (2.3.6) 
2 (¢1 + a1) To 
qu’on peut écrire en utilisant les fonctions sans dimension Ao, A}, I’, 
4 
1—Ao+ (1+ Ae) 23+ |/[0—4o)+ (1+ 4) AP — 4) 8 
: mS 


La valeur de la pulsation « de l’onde réfractée, 
pergue par un observateur qu’entraine le fluide 
dans la région EF; est alors donnée par (2.0.6) d’ou 
lon tire: 

w, 1 singo 


= : A 2.3.7 
oo I’ sing, ( ) 

Pour un observateur lié aux axes qu’entraine 
le fluide situé en Kp il y a glissement de fréquence 
et la pulsation percue sera @}. 


@,=o,{1+ =, C08 P1 
: (2.3.8) 


Y sing~o wo 
= {1+ — cosqg,) ——— —. 
( . I n) sing; I’ 
2.3.2. L’onde acoustique d’accompagnement 
eee) 
we << ¢)— a; 
Nous allons montrer que l’onde réfractée peut 


alors étre représentée par un potentiel des 
vitesses de la forme: 


Dix elias 61) 6 + lig + 0) @ + Cin’ oD ules (2.3.9) 

®, doit satisfaire 4 ’équation générale du son 
A? Dy = (1/c?) 62@1/0t? ce qui entraine les rela- 
tions: (2.3.10) 
C2924 C'2y/2— 6? —w?, 6+ 0'n'’=—01 01. 


D’autre part, sur le front S, c’est-a-dire pour 
_ x= ajt, Vexposant de ®j, doit prendre la forme 
iQ(t—y/w), dou: 


(2.3.11) 
; a /. tap hay’ t:.0% _ 2 
eter emer + 6) 10, aaa: Log 
De (2.3.10) et (2.3.11) on tire alors: 
1/2 =1/(1— Aj), (2.3.12) 


=a / au! Sep! 1(2.3.13 
pad, (2) a A?)—1, (2.3.13) 


L’éq. (2.3.13) montre que ce type d’onde n’existe 
que si la vitesse de phase w est inférieure a 


c; /1— A? = |/c2 —a?; par suite, cette onde n’ap- 
parait que dans le cas étudié au présent para- 
graphe. D’autre part, cette condition d’existence 
implique que l’on a cj > a;. Aucune onde de ce 
type ne peut donc apparaitre dans la région Eo 
ou co < ag. 

Avant de décrire onde d’accompagnement 
nous devons vérifier que la condition o <1 est 
satisfaite, c’est-a-dire que nous avons toujours 
kd/ce? <1. Remplagant dans (2.3.13) Q par sa 


valeur tirée de (2.0.6), nous obtenons: 


Ann Pn sin Yo ee . 
1 AE ec, (2.3.18) 


kqw/c? est trés petit; A; et J’ sont compris 
entre 0 et 1, on a donc kd/c? trés petit. 

Pour mettre en évidence les propriétés de 
l’onde nous prendrons des axes 0 X Y dont lori- 
gine o est entrainée par le plan S et se déplace vers 
les y croissants a la vitesse w. Nous posons donc: 


(2.3.19) 


Le potentiel ®,, prend alors la forme, indé- 
pendante de t: 


X=%x1— at, Y=y—wt. 


Dyan gt DAT ENC (2.3.20). 

L’onde est done immobile par rapport aux axes 
oX Y. Son amplitude décroit comme e& X/% Jors- 
qu’on s’éloigne de l’onde de choc (vers les x néga- 
tifs). Elle constitue donc un phénoméne qui court 
derriére l’onde de choc tout en glissant transver- 
salement a la vitesse w. Ces propriétés justifient 
le nom d’onde d’accompagnement que nous lui 
avons donné. 

Précisons sa structure. Vers les Y croissants 
le phénoméne est spatialement périodique et sa 
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longueur d’onde est: 
Ay =22¢4/n'=220/Q (2.3.21) 


ce qui était d’ailleurs évident étant donné son 
origine. Vers les X décroissants le phénoméne 
est spatialement pseudopériodique amorti. Sa 
longueur d’onde est: 


Ax =2n0e1/n = 20 (c{ —aj)/a1 Q. (2.3.22) 


Les plans d’égale phase, ou pseudo-plans d’onde 
sont donnés par: 7X+%7’ Y=const. Ils font 
done avec le plan S de l’onde de choc langle x 
donné par: 


tgy=—7'|n= (ci a;)/a,w (2.3.23) 


formule qui permet d’en donner une construction 
géométrique trés simple. 
Le taux d’amortissement de l’onde vers les X 
négatifs est: 
= 2 hw 
g c1 2 / a ee ar 


1 ci—a? 


(2.3.24) 


Il est nul pour les valeurs limites w= +/c?—a?. 
Il croit lorsque w diminue. 

Pour w=0, on a d’une part: 2=0, d’autre 
part: 

Q/w=wosingo/co et cosg~o=co/ao- 


€1/c, prend alors la valeur: 


oY Pinas? 
(:) cal al Va, Eo) 0 
aoe : : : 
C1 /w=0 Ve? me a? ao a) 
qu’on peut exprimer en fonction de la variable & 
qui définit l’amplitude de l’onde de choc, et de 


la longueur d’onde Ap = 2 2 ¢9/mo de Vonde inci- 
dente dans la région Ko: 


1) age 
£1 /w=0 Ao 


Pour calculer les caractéristiques de l’onde 
d’accompagnement a partir de celles de onde 
incidente il suffit, dans les relations (2.3.12) a 
(2.3.16), de remplacer 2 et w par leurs valeurs 
tirées de (2.0.5) et (2.0.6), savoir: 


ead 


Te: 23:25) 


2 = wo (co —ao cos yp) /co, 
w = (co — ag cos qo) /singo. 


Les formules obtenues sont compliquées. Il est 
inutile de les écrire explicitement ici. 


2.4. Constructions géométriques 


2.4.0. Construction d’Huyghens (Fig. 2) 


Le plan de la figure étant immobile par rapport 
au milieu en dépression Eo, S et S’ sont les posi- 
tions du front d’onde de choc aux époques t= 0 
et t=1. La distance SS’ est apo. 
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Considérons dans le milieu Ep, a droite de S, 
une onde plane, progressive et uniforme qui, 4 
l’époque t=0, coupe S au point Oo. A Pépoque 
t—1 cette onde coupera S’ en un point D, et son 
prolongement virtuel sera tangent a un cercle Co, 
de rayon co et de centre Oo, entiérement situé 
a gauche de S’ puisque co < ao. 

Soit B le pied de la perpendiculaire abaissée 
de Oo sur S’. BD sera égal a la vitesse de phase — 
w observée sur le front de l’onde de choc. 

Or, on peut toujours faire passer par D deux 
tangentes au cercle Co. La figure confirme donc 
quwil existe dans la région Ep deux ondes Io, Io’ 
auxquelles correspond sur S la méme vitesse de 
phase w. Ces deux ondes sont incidentes puisque 
le point de contact des tangentes menées de D a 
Co est a gauche de S: région ot Kp est virtuel. 

Considérons maintenant la région E,. Le fluide 
s’y écoule vers la droite avec la vitesse U. Une 
onde plane passant a l’époque t= 0 par le point 
Oo sera, 4 l’époque t= 1, tangente a un cercle C; 
de rayon c, et dont le centre Oj; est a la distance 
Oo One U de Oo. 

O,B est alors la vitesse a; du front d’onde de 
choc par rapport au fluide emplissant la région 
E,. Or a; <c,, par suite le cercle Cj coupe la 
droite S’ en deux points F et F’. 

Les ondes incidentes I; sont celles dont le point 
de contact avec C; est a droite de S’ (la ot Ey; 
est virtuel) les ondes réfléchies ou réfractées 
touchent C,; a gauche de S’. 

Soit maintenant D un point de S’ situé a Pex- 
térieur du segment FF’. De ce point on peut 
mener deux tangentes 4 C; dont lune correspond 
a une onde incidente Ij, l’autre 4 une onde réflé- 
chie ou réfractée Rj. 

Ainsi, si l’on se donne une onde Ry, dans Ej, 
cette onde détermine le point D, et on voit 
qu’elle peut avoir été engendrée soit par réflexion 
de l’onde I; (située dans Ej), soit par réfraction 
de lune ou l’autre des ondes Io, Ip situées dans Ep. 

Enfin, il apparait que lorsque le point D se 
trouve entre F et F’ il lui correspond encore deux 
ondes incidentes possibles Ip, Ij dans Eo, mais 
rien dans Ej; puisque D est a Vintérieur du 
cercle Cj. 

Les secteurs hachurés de la figure montrent les 
angles d’incidences go pour lesquels cette condi- 
tion est réalisée. La construction d’HuyeHENs 
ne nous permet pas alors de connaitre ce qui se 
passe dans la région Ej. 


2.4.1. Construction de l’onde d’accompagnement 
(Fig. 2) 
La formule (2.3.23) donne langle y que fait le 
pseudo-plan d’onde (plan d’égale phase) de l’onde 
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Fig. 2. 


d’accompagnement avec le plan S. Elle se préte 
a une construction simple. 

Du point F’ menons la tangente au cercle 
C;. Elle coupe la droite 0; B en G. Or, 0;G = 
0; F2/0;B = ¢;2/a;. Comme 0; B = a, on aura: 


D’autre part, D, situé entre F et F’, étant le 
point ot Vonde incidente coupe S’ a l’époque 
t=1 ona BD, =w. Par suite: 


c? /ay —ay 


7a, oa 


L’angle BD;G est donc l’angle 7 cherché. Le 
prolongement, a gauche de S’, de la droite GD 
donne le pseudo-plan d’onde de l’onde d’accom- 
pagnement, 


ONDES DE CHOC 


159 


La construction rend alors 
trés claire la variation de la di- 
rection de ce pseudo-plan d’onde 
en fonction de w, c’est-a-dire 
de la position du point Dj. 

La Fig. 2 montre ces diverses 
constructions pour toute une 
série d’incidences. Les fléches 
placées sur la trace des diffé- 
rentes ondes indiquent la direc- 
tion de leur propagation. Le 
tracé correspond au cas ow 
VYonde de choc comporte un 
doublement brusque de la pres- 
sion, c’est-a-dire ot l’on a: 


jill Sa sain 


2.5. Utilisations des phéno- 
ménes acoustiques a l’étude 
des ondes de choc 


Les conditions cinématiques 
s’expriment par des formules 
dans lesquelles n’entrent que 
des rapports de vitesse: vitesse 
du son, vitesse du front d’onde, 
vitesse d’écoulement. Elles ne 
postulent pas la validité des re- 
lations théoriques qui expriment 
ces vitesses en fonction de l’am- 
plitude du choc. 

En outre, dans les cas ot le 
front d’onde de choc est précédé 
ou suivi de phénoménes de com- 
pression ou de détente entrainant des variations 
de la température et de la vitesse d’écoulement 
du fluide, les ondes acoustiques subissent des 
réfractions continues, phénoménes par ailleurs 
bien connus. 

Par suite, l’étude expérimentale de la ciné- 
matique des ondes acoustiques réfléchies ou ré- 
fractées par un front d’onde de choc fournit une 
méthode de recherche qui n’est probablement 
pas sans intérét pour les aérodynamiciens. On 
peut envisager pour la mettre en ceuvre, soit des 
procédés photographiques, voire cinématographi- 
ques montrant la forme et les déformations des 
fronts d’onde, soit des méthodes microphoniques 
décelant les variations de fréquence. 

On remarquera que nombre d’informations 
peuvent étre recueillies par la seule étude de la 
réflexion des ondes, c’est-a-dire par des mesures 
effectuées uniquement dans la région en com- 
pression Ej. 
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3. CONDITIONS DYNAMIQUES 


3.0. Elimination des ondes thermiques 


Partons des relations établies au paragraphe 
1.1.2, savoir: 


6; -—09 =2€& = (@1 — @o), (11.17) 
ity a 
S$, — So = é aha (@1 — @0)s (1.1.18) 


i oe 
Ci 1 il 


¢1 uy —G a 


Nous avons vu que, dans notre probléme, la 
région en dépression Eo ne peut contenir que des 
ondes incidentes, représentables par le potentiel 
des vitesses: 


(« Peay: ex} (1.1.21) 
y 


nes eivelt— (£9 COSPo + Yo SiN Po)/Co} (3.0.1) 
On aura donc: 
y= —i(wo/co?) Po, 
Wy =e c?) Do, 
Oo = Y So iy (@0/¢) Po (3.0.2) 
ug = —i(o9]¢0) eos 70 Bo = co cos 70-5. 


Dans la région Ej, il y a superposition d’un 
phénoméne acoustique (lettres affectées de l’in- 
dice 1a) et d’un phénoméne thermique (indice It). 

Les formules (2.1.7) nous donnent pour le 
phénoméne thermique 


@14=0, O4¢=—Sit, Wie = —i(kQ/y ay) 514, (3.0.3) 
Dans le phénoméne acoustique nous avons: 


(3.0.4) 


O14 = Y Sla 


et au total: ©} = @j;4+ O1g = @1q, 01 = 614+ O14; 
$1 = Sit + Sia, U1 = Ue + Uig- 
L’équation (1.1.18) donne alors: 


Sig = —y J (@1 — Wo) (3.0.5) 


ou J est la fonction définie par la formule (1.1.11). 
Cette relation nous donnera l’onde thermique 
lorsque nous aurons calculé l’onde acoustique. 
Portant cette valeur de sj; dans la troisiéme 
relation (3.0.3), nous avons: 


u1p=i(kQ/a1)J(@1—@o) _—(3.0.6) 


que nous portons dans (1.1.21), qui devient: 


Cl Ujg (c 1 


aick ; ae ye I 
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Or, nous savons que les nombres k2/ajc; et 
k Q/a, co sont trés petits. D’autre part, le Tableau 
II du Chapitre 1 nous montre que la fonction J 
est toujours inférieure aux fonctions G et H. Par 
suite, dans l’équation (3.0.7) nous pouvons négli- 
ger les termes imaginaires, c’est-a-dire les termes 
provenant de l’onde thermique. 

Il nous reste donc une condition qui n’intéresse 
que les ondes acoustiques existant de part et 
d’autre du front S de l’onde de choc. 

e @1 ta @o0 


Y 


). (3.0.8) 


Pare G =F (co wo —H 


3.1 L’onde réfléchie 


Nous avons a considérer dans la région EK; 
deux ondes, l’une incidente, l’autre réfléchie, re- 
présentées par leurs potentiels des vitesses, 


(3.1.1) 
(3.1.2) 


D’autre part, la région Eo ne contient aucune 
onde. Le second membre de la condition (3.0.8) 
est donc nul. 

Supposant les conditions cinématiques satis- 
faites (3.0.8), donne alors: 


B, 1 (cosy; — G) + Bi w1( cosg; — G) = 90. (3.1.3) 


iw, [t— (x, cosg, + ¥: Sin g,)/¢e) 
@®,=— Bye 1 1 1 1 1 oe 


PD 2M Be elo [t— ,coser + ys sin gy’ )/es) 


L’amplitude des ondes est proportionnelle aux 
facteurs Bw. Le rapport de l’amplitude de onde 
réfléchie a l’amplitude de l’onde incidente, c’est-a- 
dire le coefficient R de réflexion des amplitudes, 
est alors: 


R = —(cosqy1 — G)/(cosy, —G). (3.1.4) 


Reportons-nous a la relation cinématique 
(2.2.1) qui lie les tangentes 11, 1; des demi- 


1 1 Pe: 
angles 5 71, 5 i, nous pouvons éliminer yj, et 


il vient: 
1 AN 1 Cae (3.1.5) 
+73 =~ 
NEE AS 1+6G 
l+c ee AY, 2 jr ei 
14+ Ay 1+6 I 


avec tj7<(1—Aj)/(1+ Aj) (onde incidente). 

Pour gy; = 0, t; = 0 (incidence normale) Ro 
=(1—G)/(1+ G), négatif car G>1 (voir Ta- 
bleau II, Chapitre 1). 

Pour ¢ = arccos Aj, tT) = +V(1—A))/(1+ Ay) 
(incidence limite). Rjim, = —1. 

Comme, en ce cas, on a 91 = G1, | = 1, 
cette valeur signifie que, dans la région ot les 
deux ondes se superposent (et qui, d’ailleurs, est 
évanouissante) leur somme est nulle. En d’autres 
termes, le phénoméne ne comporte qu’une onde 
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acoustique incidente qui s’arréte sur S ce qui est 
bien le cas. 

Le coefficient de réflexion des amplitudes est 
done constamment négatif, et, en valeur absolue, 
eroit de (G—1)/(G+ 1) a1 quand ¢q, croit de 0a 
arccos A. 

L’onde réfléchie ne sera donc pratiquement 
observable sous des incidences voisines de la 
normale que si l’onde de choc est intense (G 
_nettement supérieur a 1). 


3.2. L’onde réfractée progressive 


Représentons l’onde incidente, située dans Eo, 
par son potentiel des vitesses 


Do aes Bo a @o [E — (Xo COS Go + Yo SiN Yo)/Co] 


_et l’onde réfractée, située dans E), par 


= Bi elo lt— (®, COSp: + ¥1 Sin y:)/¢1) 


Dés Vinstant que les conditions cinématiques 
sont satisfaites les exposants des exponentielles 
sont les mémes sur S. La condition dynamique 


(3.0.8) donne alors: 


IB, @1(cosgy;— G) = Bowo(cosgo—H). (3.2.1) 


Nous comparerons les amplitudes des 
pressions acoustiques dans les deux ondes, 
e’est-a-dire que nous formerons le rapport: 


R = P10 1/ Po Ho = 01 Bi o1/00 Bo wo 
dont (3.2.1) donne la valeur: 


Fa 1 H—cosgo 
~ RI G—cosqi 
ot R= go/o1 et 
phd sib e+e af 
RP} €(1+ p46) é 
(voir Tableau II, Chapitre 1). 
Pour calculer RF il faudra se reporter aux condi- 
tions cinématiques donnant qj, en fonction de qo. 
Nous donnerons seulement ci-dessous les 
valeurs R en fonction de £ pour go = 0 (onde 
_incidente paralléle A l’onde de choc et fuyant 


devant cette derniére) et go = 2 (onde incidente 
paralléle 4 onde de choc et venant 4a sa ren- 


(3.2.2) 


contre). 
Nous avons dans les deux cas yj =2, d’ou: 
S71 
eeelaee Romar CTT 
1 H+1 
NSE ee eo 


RP el 
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1 0 1 
0,9 0,000 2 1,094 
0,8 0,0010 | pay 
0,7 0,0035 1,36 
0,6 0,009 2 55 
0,5 0,022 1,82 
0,4 0,049 2,19 
0,3 0,120 2,79 
0,2 0,332 3,92 
0,1 1,29 7,00 
0,05 3,85 LISA 
0,025 7,85 17,9 
0 fore) fone) 


3.3. L’onde d’accompagnement 


L’onde incidente étant représentée comme pré- 
cédemment, l’onde d’accompagnement dérive du 
potentiel 

QD; = By elias —06i:)t+ [(in + 6) a. + inf'yller 


Lorsque les cenditions cinématiques sont satis- 
faites, la condition aux limites (3.0.8) donne: 


(3.3.1) 
: H —cosqo ¢ in+€ 
iw@o By ————— _ = (iwi —4}) Bi | G—- —— 
co i@,— 01 


Le rapport R qui nous intéresse le plus est 
celui de la pression acoustique immédiatement 
derriére onde de choc py = —¢@1(im1—4}) By, a 
la pression acoustique dans l’onde incidente 
po= —1 00 0 Bo . 

Tenant compte des relations (2.3.12 a 2.3.17) 
on obtient: (3.3.2) 


cos~o — H i—0,/@ 


~ LPR i(4;—G6)— (6/1) 1/41 — 6) 


et, du fait que 6;/@;= A, Video A) (cy/w)2—1: 


ey Wee era Cie ethan 
cet igh SH (I— 4, 6)2+(G2—1)W2 


ou l’on a posé: 


a ae et NESE 


Cy sin 0 


A titre d’indication nous avons effectué le cal- 
cul de |R| dans le cas extréme o w= 0, c’est-a- 
dire cosgg= 1/Ag. On a alors: 

RI Pen TAs 
sor piritecen egy Pe RRA 


Pour go = arecos1/Ag on obtient 


£=10,9 0,8 0,7 0,6 0,5 0,4 0,3 0,2 0,1 0,05 0,025 0 
|R| = 11,05 1,15 1,31 1,46 1,71 2,06 2,66 3,86 7,45 14,6 22,8 oo 
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On voit que, comme dans le cas précédent, les 
pressions acoustiques peuvent étre majorées dans 
une proportion importante. 


4. RESUME DES RESULTATS OBTENUS 


La présence d’ondes acoustiques dans l’une ou 
Vautre des régions que sépare le front d’onde de 
choc se traduit par les phénoménes suivants: 

a) Dans la région en dépression Eo, dans la- 
quelle onde de choc progresse a une vitesse 
supersonique, aucune sorte d’onde réfléchie ou 
réfractée ne se produit. 

b) Dans la région en surpression Kj, la ren- 
contre d’une onde acoustique incidente (pro- 
venant de la région Eo ou de la région Ej) et de 
Yonde de choc provoque la formation de deux 
ondes de type différent. L’une est une onde 
acoustique (réfractée ou réfléchie selon l’origine 
de l’onde incidente), l’autre est une onde ther- 
mique. 

c) L’onde thermique, en premiére approxima- 
tion (c’est-a-dire lorsque les ondes acoustiques 
incidentes sont d’une amplitude usuelle), n’est 
pas progressive. Elle est done immobile par 
rapport au gaz contenu dans la région en sur- 
pression et entrainée par ce gaz dans son mouve- 
ment d’écoulement d’ensemble. Elle s’éloigne 
alors du front de l’onde de choc. Cette onde com- 
porte des variations périodiques spatiales de tem- 
pérature et de densité, mais aucune variation de 
pression. Elle s’amortit lentement et apériodique- 
ment. 

d) Lorsque Vonde acoustique incidente est 
située dans la région en surpression, il y a produc- 
tion dans cette méme région d’une onde acousti- 
que réfléchie du méme type que l’onde incidente, 
c’est-a-dire dans le cas étudié ici, plane, pro- 
gressive et uniforme. 

La construction d’Huyc¢HENs permet d’en ob- 
tenir les caractéristiques cinématiques: direction 
de propagation, longueur d’onde. 

L’amplitude des variations de pression dans 
Vonde réfléchie est toujours inférieure 4 l’ampli- 
tude de ces variations dans l’onde incidente. 

e) Lorsque Vonde acoustique incidente est 
située dans la région en dépression il y a produc- 
tion, dans la région en surpression, d’une onde 
acoustique réfractée dont les propriétés dépen- 
dent de la vitesse de phase w du phénoméne 
engendré sur le front du choc par l’onde inci- 
dente. Les propriétés principales de l’onde réfrac- 
tée sont les suivantes: 

x) Il existe toujours dans la région en dépression 
Eo deux incidences d’ondes planes qui pro- 
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duisent sur le front de choc la méme vitesse 
de phase w, et par suite, dans la région en 
surpression E; la méme onde réfractée (au 
rapport d’amplitude prés). 

L’observation de lVonde réfractée dans la 

région en surpression E; ne permet done pas 
de discerner si cette onde a été engendrée par 
lune ou l’autre, ou toute combinaison, des 
deux ondes incidentes possibles. 
Si la vitesse de phase w de la perturbation 
engendrée sur l’onde de choc par londe 
acoustique incidente est supérieure a une 
certaine valeur limite wo (qui dépend de 
Vintensité du choc), onde réfractée est du 
méme type que l’onde incidente: plane, pro- 
gressive et uniforme. La construction d’Huy- 
GHENS permet d’en obtenir les caractéristi- 
ques cinématiques. 

L’amplitude des variations de pression dans 
Vonde réfractée peut étre selon les cas (angle 
d’incidence, amplitude du choc) inférieure, 
égale ou supérieure a l’amplitude de ces varia- 
tions dans l’onde incidente. Toutes choses 
égales d’ailleurs, cette amplitude croit avec 
amplitude du choc. 

La majoration d’amplitude peut devenir 
trés importante dans le cas de choes intenses. 
Si la vitesse de phase w est inférieure a la 
vitesse limite wo, ce qui correspond a des 
incidences situées dans un secteur dont les 
limites sont fonction de amplitude du choc, 
Vonde réfractée prend la forme d’une «onde 
d’accompagnement». 

Cette onde reste attachée a l’onde de choe 
qu’elle accompagne tout en glissant paralléle- 
ment a cette onde a la vitesse w. 

Les plans d’onde (plans d’égale phase mais 
non d’égale intensité) sont obliques par rap- 
port au choc. Leur direction peut étre ob- 
tenue par une construction géométrique 
simple. 

L’amplitude de Vonde décroit exponen- 
tiellement lorsqu’on s’éloigne du front de 
Ponde de choc. Immédiatement derriére 
Ponde de choc, amplitude des variations de 
pression dans londe d’accompagnement est 
d’autant plus grande que le choc est plus 
intense. Elle peut étre supérieure a l’ampli- 
tude des variations de pression dans londe 
incidente. 

En résumé l’existence d’ondes acoustiques 
dans la région en dépression, dans laquelle 
progresse onde de choc, peut se traduire 
derriére Ponde de choc d’une part par un 
champ d’ondes acoustiques progressives qui 
se répandent peu a peu dans la région en sur- 
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pression toute entiére, d’autre part par l’appa- 
rition d’ondes d’accompagnement, c’est-a- 
dire de perturbations oscillantes que l’onde 
de choc entraine derriére elle dans sa pro- 
gression. 

Les unes et les autres peuvent dans certains 
cas comporter des oscillations de pression trés 
supérieures a celles qui existaient avant 
Varrivée de l’onde de choc. 

(Regu le 11 Décembre 1954.) 
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ANZIEHENDE UND ABSTOSSENDE KRAFTE ZWISCHEN KUGELN 


IM SCHALLFELD 
von W. DORR 


III, Physikalisches Institut der Universitat Géttingen 


Zusammenfassung 


Es wurden die Kriifte, die in einer stehenden Schallwelle zwischen zwei Kugeln auftreten, 
gemessen. Und zwar fiir Abstinde, die mit den Durchmessern der Kugeln vergleichbar sind, bei 
denen also die Kénr1csche Formel (1) nicht mehr giiltig ist. Aus den Messungen ergab sich ein 
Zusatzglied AK zur K6nicschen Formel (3), das in seinem Aufbau mit dem theoretisch geforder- 
ten (2) iibereinstimmt, soweit es sich um die Abhingigkeit vom Abstand, von der KugelgréBe 
und der Schnelle handelt. Nicht iibereinstimmend ist die Form der Winkelabhangigkeit: experi- 
mentell ergibt sich eine Zunahme von AK mit 0 (Winkel zwischen Kugelverbindungslinie und 
Stromrichtung), wihrend theoretisch 4K mit wachsendem § abnehmen miifte. 


Summary 

The forces that arise between two spheres in a stationary sound wave are measured, especially 
for distances between the spheres comparable with their diameters, for which K6nic’s formula 
(1) ist not valid. From the measurements, a correction 4K is found for K6nr1c’s formula (3) which 
agrees in form with a theoretically derived formula (2), in respect of dependence on distance, 
sphere size and velocity amplitude. The variation with the angle 0 between the line joining the 
spheres and the axis of motion is not in agreement; the experiments would make AK increase, 
the theory decrease with 6. 


Sommaire 

On a mesuré les forces prenant naissance entre deux sphéres, dans une onde sonore stationnaire, 
dans le cas ow les distances sont comparables aux diamétres des sphéres, et ou, par conséquent, la 
formule de K6nic (1) n’est plus valable. Il résulte des mesures qu’il faut ajouter a la formule 
de KOni¢ un terme additionnel 4K (3), dont la constitution est conforme 4 celle qu’exige la 
théorie (2), du moins en ce qui concerne l’influence de la distance, de la grosseur des sphéres et de 
la vitesse. L’accord n’est pas satisfaisant, par contre, en ce qui concerne l’influence de l’angle: 
d’aprés l’expérience, 4K augmente avec @ (angle que fait la droite joignant les sphéres et la 
direction du courant), tandis que, d’aprés la théorie, AK devrait diminuer quand @ augmente. 


Einleitung 

Befinden sich feste Kérper in einem bewegten 
Medium, so treten zwischen ihnen hydrodynami- 
sche Krafte auf. Diese Krafte spielen wahr- 
scheinlich eine Rolle bei der Entstehung der 
Kundtschen Staubfiguren (Rippelbildung) sowie 
bei der Koagulation von Aerosolen durch starke 
Beschallung. Ihre GréBe theoretisch anzugeben, 
wird im allgemeinen Fall schwierig sein. Nur bei 
einfachen Kérpern (Kugel) bleibt die Rechnung 
iibersichtlich, fiihrt allerdings auch nur zu mehr 


oder weniger guten Naherungen. Denn einerseits 
ist in der hydrodynamischen Rechnung das 
Geschwindigkeitspotential der Strémung um 
zwei Kugeln nur naherungsweise gegeben, etwa 
in Form einer Entwicklung, deren héhere Glieder 
vernachlassigt werden; wahrend andererseits der 
Einflu8 der Reibung véllig auBer acht gelassen 
wird. Das ist in Luft durchaus gerechtfertigt, so- 
lange der Oberflachenabstand der Kugeln gréBer 
ist als die Dicke der Zahigkeitsgrenzschichten an 
den Oberflachen. Wird der Abstand jedoch kleiner, 
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so wird man das nicht mehr tun diirfen, sondern 
man muf die durch Zahigkeit bedingten Effekte 
mit in Rechnung stellen. 

Eine Formel fiir die Krafte, die in einem 


Schallfeld zwischen zwei Kugeln wirken, hat’ 


Kénie [1] angegeben. Man kann sie in folgender 
Form schreiben: 
3 rene 


3 1 
Io y UL 5a a (5 cos2 6+ 3) 


(1) 


Taiis : 
K,= 370; i v2 sin 8 cos. 


Darin bedeuten: K; = Kraft in Richtung der 
Kugelverbindungslinie, K,—= Drehmoment des 
Kugelpaares, rj, r2 = Kugelradien, 6 = Zentral- 
abstand der Kugeln, 6 = Winkel 
Kugelverbindungslinie und Stromrichtung, v9 = 
Spitzenwert der Schnelle v=—vosinwt, 0, = 
Dichte der Luft. 

Der Rechnung liegt die Annahme eines inkom- 
pressiblen und reibungslosen Mediums zugrunde; 
auBerdem ist die Stérung, die der einfallende 
Strom am Ort der ersten Kugel durch den von 
der zweiten Kugel ausgehenden Reaktionsstrom 
erfahrt, vernachlassigt. Die Formel gilt daher 
nur in dem Bereich, in dem obige Vernach- 
lassigung zulassig ist, also dann, wenn der Ab- 
stand 6 >1r,+ 19 ist. 

Bei kleineren Abstanden muf man diese Sté- 
rungen héherer Ordnung mit bericksichtigen. 
BJERKNES [2] hat einen solchen Fall berechnet. 
Danach ergibt sich ein Zusatzglied AK zum 
Ko6niéschen Wert K,. Es ist dann die Gesamt- 
kraft: 


zwischen 


mit (2) 


3 3 3 1 
AK=—~ xo, rr net a7 (1 + 3 cos? 6) 05. 


2 


Diese Formel gilt auch fiir Abstande 6 <r; + re; 
allerdings mit der Einschrankung, da man eine 
richtige Beschreibung nicht mehr erwarten kann, 
wenn Effekte, die durch die Zahigkeit des 
Mediums bedingt sind, auftreten. 


Messungen 


Von den SchallfeldgréBen geht nur die Schnelle 
v in die Formeln (1) und (2) ein. Die Messungen 
erfolgten daher in einem Kundtschen Rohr im 
Schnellebauch. Und zwar wurde zur Erreichung 
hoher Schallstarken eine von OBerst [3] ange- 
gebene Anordnung benutzt. Sie besteht aus 
einem Rohr mit sprungartiger Querschnitts- 
verengung, das am Anfang des weiteren Teiles 
erregt wird. An der Sprungstelle tritt dann eine 
Druckiiberhéhung ein, die wiederum als An- 
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fangsdruck auf den engeren Teil wirkt und so in 
ihm héhere Schallstarken hervorzurufen vermag. 
Die Erregung erfolgte durch einen 10 Watt-Laut- 
sprecher, betrieben tiber einen Kraftverstarker 
von einem Philips-Schwebungssummer, 

Gemessen wurde bei f = 525 Hz die rein an- 
ziehende oder abstoBende Kraft K, zwischen 
Kugeln verschiedener GréBe (r = 0,44; 0,37; 
0,28 cm). Die Kugeln waren aus diinnwandigem 
Glas geblasen und hingen je an zwei gespreizten 
Faden in den Schnellebauch des kleinen Rohres 
hinein. So wurde die Wirkung des Drehmomentes 
Ke ausgeschaltet. Die Kraft ergab sich aus der 
mikroskopisch gemessenen Auslenkung der Ku- 
geln gegen ihre Ruhelage. Die so erreichte stati- 
sche Empfindlichkeit ging bis zu 1,25 dyn/em. 
Bei mikroskopischer Ablesung sind also Krafte 
bis herab zu 1,5-10-3 dyn noch mefSbar. 

Die Schnelle im Schwingungsbauch bestimmte 
sich aus dem Schalldruck am Ende des Rohres, 
was durch Amplitudenmessungen von Schwebe- 
teilchen direkt am Ort der Kugeln kontrolliert 
wurde. 

Effekte, welche die zu messenden Anziehungs- 
krafte stérend tberlagern kénnten, waren nach 
Méglichkeit auszuschalten. Als solche kommen 
in Frage: 

1. Wandeinfliisse, 2. Schallstrahlungsdruck, 
3. Turbulenz, 4. Stokessche Reibung. 


1. Die Wandeinflisse, soweit sie hydro- 
dynamischer Art sind oder von der akustischen 
Grenzschicht herriihren, lassen sich unwirksam 
machen, wenn der Rohrquerschnitt grof genug 
gegen die Abmessungen der Kugeln ist. Das war 
bei einem Querschnitt von 5 * 5 cm? der Fall. 
Die Zirkulationsstrémung .im Rohr schlieBlich 
kann klein genug gehalten werden, wenn man 
mit der Schnelle nicht zu hoch geht. Bei Schnelle- 
werten von 40 cm/s trat sie noch nicht stérend in 
Erscheinung. 

2. Der Schallstrahlungsdruck trat nicht 
auf, da im Schwingungsbauch gemessen wurde. 

3. Turbulenz. Die Formeln iiber die hydro- 
dynamischen Krafte zwischen Kugeln gelten 
natiirlich nur so lange, wie die Strémung um die 
Kugeln laminar ist. Nach MOLLER [4] ist das — 
auch noch im Gebiet hinter der Kugel — der Fall, 
solange die Reynoldszahl Re < 300 ist. 

Es gilt also fiir eine Kugel vom Durchmesser 
d=1 cm in Luft: 


d.h. es muB v < 43 cm/s sein. Damit ist eine 
obere Grenze fiir die Schnelle gesetzt. 
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4. Durch die Stokessche Reibung greift 
eine Kraft an den Kugeln an, die sie veranlaBbt, 
an der Bewegung des umgebenden Mediums zu 
einem gewissen Grade teilzunehmen, wodurch 
sich die Relativschnelle Kugel—Luft verandert 
(Mitnahme). Den Grad der Mitnahme in vor- 
liegendem Falle theoretisch exakt anzugeben, 
wire einigermaBen kompliziert und zwange zu 
umfangreichen Betrachtungen. Denn die Messun- 
gen spielen sich bei Reynoldszahlen von 100---300 
ab, wahrend die einfache Stokessche Formel bis 
Re = 0,4 giltig ist und auch die Oseensche 
Erweiterung nur bis zu Re = 2 reicht. 

Nun sind aber die Kugeln bei der Messung 
nicht frei beweglich, sondern mit der Eigenfre- 
quenz f= 1 Hz pendelnd aufgehangt. Die Schall- 
frequenz betragt etwa 500 Hz. Sollen die Kugeln 
mitschwingen, so miibten sie erzwungene Schwin- 
gungen von 500 Hz ausfiihren. Um dabei aber zu 
wirksamen Amplituden zu gelangen, mite die 
anregende Kraft von einer GréSenordnung sein, 
die im ganzen System nicht auftreten kann, wie 
eine einfache Abschatzung zeigt. Im tibrigen kann 
man recht gut mikroskopisch kontrollieren, ob 
eine Mitfiihrung der Kugeln stattfindet. Bei 40- 
facher VergréBerung und Luftamplituden von 
0,28 mm wurde keine meBbare Mitfiihrung fest- 
gestellt. 

Es folgen nun einige MeBergebnisse: 

Halt man 6 konstant, so ergibt sich nach Gl. (1) 
ar (2: K; = const - v3 

Andert sich 6 wahrend. der Messung nur wenig, 
so gilt noch: as 
Abb. 1 zeigt eine Messung dieser Art. Bei 0 = 0 
wurde die Schnelle von vp = 0 bis vp = 35 cm/s 
gesteigert und die zu den jeweiligen Schnelle- 


1000 cm?/s? 


0 200 400 #600 800 


i 


Abb. 1. Abhangigkeit der Anziehungskraft von der Schnelle 
6 = 0, Kugeln: r = 0,74 em. 
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werten gehérigen Abstinde der Kugeln wurden 
gemessen. Der dabei iiberstrichene Abstands- 
bereich war sehr gering (2% des Gesamtabstan- 
des), so das nach den Formeln (1) und (2) die 
Beziehung K,64= const-v? zu erwarten war. 

Wie man sieht, wird dies experimentell recht 
gut bestatigt. K, hangt also quadratisch von der 
Schnelle ab. 

Weiterhin wurde die Abhangigkeit der Kraft 
vom Abstand untersucht. Bei konstanter Schnelle 
und festgehaltenem Winkel 0 wurde 6 verandert 
und die zu jedem Abstand gehérige Kraft ge- 
messen. Abb. 2 zeigt diesen Zusammenhang. Ky 
ist die gemessene Kraft, K, die nach Formel (1) 
berechnete. 

Es ergibt sich, daB mit kleineren Abstanden 
die Differenz K,— Ky immer gréer wird, wie 
es ja auch nach den Voraussetzungen der 
Formel (1) zu erwarten ist. 


0,12 
dyn 


0,08 


0,04 


{7 08 09 


ey) 
Abb. 2. Abhangigkeit der Anziehungskraft vom Abstand 0; 
6 = r/2, Kugeln: 7, = 0,375 cm, rz = 0,345 cm, 
Vv. = 38 cm/s. 


40cm 


SchlieBlich wurde auch noch die Winkelab- 
hangigkeit der Kraft gemessen. Bei konstanter 
Schnelle und kleinen Abstanden wurde die Ver- 
bindungsachse der Kugeln um 7/2 gedreht und 
von 10° zu 10° die Auslenkung der Kugeln gegen- 
iiber ihrer Ruhelage bestimmt. 

Abb. 3 zeigt als Ergebnis einer solchen Messung 
Ra Br(9) 96 

K, (0) 04 
vom Winkel 6 aufgetragen. Die gestrichelte Kurve 
gibt den Wert von R, wie er aus der KOn1Gschen 
Formel (1) folgen wiirde. 

Wie zu erwarten, tritt auch hier eine Abwei- 
chung auf, denn die Abstande sind so klein, daf 
der Giiltigkeitsbereich der Gleichung (1) iiber- 
schritten wird. 

Man kann nun aus Messungen der Art von 
Abb. 2 und 3 ein Zusatzglied AK so bestimmen, 
dafS K,—K,+ AK die gemessenen Krafte be- 
schreibt. 


das Verhaltnis in Abhangigkeit 
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"0 30° 60° 90° 


—- 6 


Abb. 3. Abhangigkeit der Anziehungskraft vom Winkel 6, 
Kugeln: r = 0,74 em. 


Anhand von Messungen, die mit drei verschie- 
denen KugelgréBen ausgefiihrt wurden, ergab 


sich, da dieses Zusatzglied fir alle Winkel 0 
negativ (anziehend) ist und sich darstellen laBt 
durch: 

2000 


1000 


06 08 10 12 


Abb. 4. x x ry = 0,485 cm, r2= 0,495 cm, g=0 u. 
AA %=0,375 em, rz= 0,345 cm, 0=0 u. x/2 
©0 T= 0)275'em, 7 —05295'em., 6 — 0 en 
Vy = 38 cm/s. 
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Sas 
Ty +13 


2 


. 74 sin? 6 + 6) v3. 
(3) 


3 
AK =— yz moL ima 


Das zeigt Abb. 4. Hier ist 


10 AK 
eo =a 2 
AK* = 378 + ey (Asin? 0) 


fir verschiedene KugelgréBen und verschiedene 
Winkel in Abhangigkeit von 6 aufgetragen. Die 
Punkte sind MeBpunkte, die ausgezogene Kurve 
ist nach Gl. (3) berechnet. 

Vergleicht man das experimentell gewonnene 
Zusatzglied (3) mit dem theoretisch berechneten 
(2), so zeigt sich, daB beide Ausdriicke iiberein- 
stimmen, soweit es die Abhangigkeit vom Ab- 
stand 6 und den Kugelradien r betrifft. Dagegen 
unterscheiden sie sich in der Form der Winkel- 


abhangigkeit : 


2 
theoretisch: AK = const <t (1 + 3 cos? 6), 


2 
experimentell: 4K = const - (4 sin? 6 + 6). 
Theoretisch miBte also AK mit zunehmendem 
Winkel 0 kleiner werden, experimentell jedoch 
ist unzweifelhaft ein Anwachsen von AK mit 6 
festzustellen. 

Herrn Prof. E. Meyer danke ich fiir die An- 
regung zu dieser Arbeit sowie fiir wertvolle Hin- 
weise wahrend ihrer Durchfiihrung. 

(Eingegangen am 16. Dezember 1954.) 
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ACOUSTIC RADIATION PRESSURE ON A COMPRESSIBLE SPHERE 


by K. YostoKa and Y. KAWASIMA 


Institute of Scientific and Industrial Research, Osaka University 


Summary 


Calculations are given of the acoustic radiation pressure on compressible spheres suspended 
freely in a plane progressive or plane stationary sound wave field in a non-viscous fluid. With 
plane progressive waves the radiation pressure on bubbles shows a prominent peak corresponding 
to resonance in the pulsation mode. With plane stationary waves, on the contrary, the radiation 
pressure on the resonant bubble vanishes and the directions of the pressure are such that the 
bubbles smaller than the resonant are forced towards the nodes and the larger ones towards the 
loops. On small spheres, of which the densities are of the same order as that of the surrounding 
medium, such phenomena in the radiation pressure disappear and the formulae contain the 
density-compressibility factors instead of Kinc’s density factors. 


Sommaire 


On calcule la pression de radiation acoustique qui s’exerce sur des sphéres compressibles sus- 
pendues librement dans un champ d’ondes sonores planes, progressives ou stationnaires, produit 
dans un fluide non visqueux. Dans le cas des ondes planes progressives, la pression de radiation 
sur les bulles présente un maximum accusé, qui correspond 4 une résonance dans le mode pulsa- 
toire. Dans le cas d’ondes planes stationnaires, au contraire, la pression de radiation s’annule sur 
les bulles en résonance, et les directions de cette pression sont telles que les bulles plus petites que 
les bulles en résonance sont déplacées vers les nceuds, tandis que les bulles plus grandes sont 
déplacées vers les ventres. Pour des petites sphéres dont la densité est du méme ordre de grandeur 
que la densité du milieu ambiant, il n’y a plus de phénoménes de ce genre, et les formules renfer- 
ment les facteurs densité-compressibilité au lieu des facteurs densité de Kine. 


Zusammenfassung 

Der Schallstrahlungsdruck auf kompressible Kugeln in einer nichtviskosen Fliissigkeit wird fiir 
fortschreitende und stehende ebene Wellen berechnet. Bei fortschreitenden Wellen ergibt sich bei 
der radialsymmetrischen Eigenresonanz ein Spitzenwert fiir den Strahlungsdruck, wahrend er in 
stehenden Wellen bei der Resonanz ganz verschwindet und kleinere Blasen sich in den Schwin- 
gungsknoten, gréBere sich in den Bauchen sammeln. Im Falle sehr kleiner Blasen, bei denen die 
Dichte gré6BenordnungsmaBig gleich der Dichte des umgebenden Mediums wird, treten diese Er- 
scheinungen nicht mehr auf, und die Formeln enthalten eine Dichte-Kompressibilitats-Beziehung 


anstatt des Dichtefaktors nach K1INc. 


1. Introduction 


The acoustic radiation pressure on rigid spheres 
freely suspended in a non-viscous fluid was cal- 
culated by L. V. Kine [1]. The result of the cal- 
culation stands in good agreement with the later 
observations [2] made in liquid on balls of such 
a hard material as steel. 

Disagreement may, however, be expected for 
spheres of softer materials, of which the com- 
pressibility effect is no longer negligible. The 
experiments by K. SOLLNER and C. Bonpy [3] on 
emulsions showed that globules less dense than 
the medium of dispersion were forced to accumu- 
late at the nodes and denser ones at the loops in 
a plane stationary wave field, the critical density 
ratio predicted by K1Nc’s theory being 0.40. The 
magnitude of the disagreement is serious for 
bubbles, on which the radiation pressure was 
recently observed by the present authors [4] to 


be 10°.--108 times as large as that given by 
Kine’s theory. 

In the following, calculations are given of the 
acoustic radiation pressure on spheres of finite 
compressibility placed freely in a plane pro- 
gressive or plane stationary sound wave field in 
a non-viscous fluid. 


2. The radiation pressure on a moving boundary 


The pressure variation in a perfect fluid, where 
the train of sound waves is propagating, is given 
to terms of the order uj/c? by 


; a ee Cae 
sp = 00 (b, + ,)— Font+ 5 26%, (1) 
where 9 is the equilibrium density of the 
medium, c the sound velocity, u; the first order 
particle velocity (= —\/®;), ® the first order 
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velocity potential given by the approximate wave 
equation /\ ®, = (1/c?) 6?@,/ét? and @2 the 
second order velocity potential. 

Let a closed boundary S in the medium be in 


motion with small velocity of the same order as - 


u; and S(t) be its position at the time t. The 
acoustic radiation pressure on the moving bound- 
ary <P» is given (to terms of the second order) 


by 
Py = | fopn df> (2) 


S(t) 


= fe et aie te 


s(t 
+<f freminap ef fag 


So indicating the boundary at its equilibrium 
position. 

If a stationary state is established throughout 
the first and second order acoustic fields, the 


of Pinal. 


second order particle velocity ug and 09 Py» will 
take the following form: 


ug = —\/ ®2 (4) 
= const + fy + Y Af. ncosnat + fo, nsinnot}, 
; i| 1 A 

20 D2 = dp2 9 20 uy 9 o Di (5) 


= const + go + )'(g1,ncosnwt+ go nsinnwt}, 
(pial! 


where dp2 is the second order pressure variation, 


fo: fi,n f2, n» 80> §1,n» §2,n functions of the co- 
ordinates only. 
From (4) and (5), we have 


8o(x, y, 2) = const, (6) 


from which follows 
<f [eoP2ndf>=0. (7) 
Sy 


Expressing the normal velocity of the boundary 


by uy, n, we have 
dt i {/ feotim ay} — J [oor maf + 


S(t) S(t) (8) 
+f foov% Un Of. 
S(t) 


If both ®, and the displacement of the boundary 
are periodic with regard to time, it follows, by 
taking the time average of (8) and writing So in 
place of S(t) in the last integral, that 


K. YOSIOKA and Y. KAWASIMA: RADIATION PRESSURE ON SPHERE 


ACUSTICA 
Vol. 5 (1955) 


<f food: mdfy = Bae 9 UE 
S(t) 


= Ulla ec un df>, 
So 


u;t being the tangential particle velocity in the 
medium. 


Putting (9) and (7) into (3), we have finally 


ame (un a + ut t) Un df> + (10) 


“ffewnny—ff3' ce ° Din df». 


The radiation pressure on any moving bound- 
ary is thus calculated in terms of the first order 
velocity potential, all the integrals being taken 
over the surface at its equilibrium position. 

When the boundary is spherical (radius a) and 
the acoustic field axially symmetrical with regard 
to the direction 9=0, the radiation pressure 
along = 0 is given by 


<P> = ¢Pw> + <Puw> + < Pur, wo + <Po>, 


where 


& 2 
(Pi = <— 7 a oof (3) sin 9 cos 8 d@>, (12) 
0 =a 


(11) 


/ Fay I P OP) ; i 
Pega cof ( a0 ). Sint coe 
oD, 0D, 
Gein uo) = Les aaaioof (Ge oy Be (Se). 
x sin?6 d@>, (14) 


(Pp = c= 


ma200 (,< : 
“| (P1)3-asin 8 cosOd0>, (15) 

0 

®, being the solution, satisfying the boundary 


conditions at its equilibrium position, of the 
scattering by a compressible sphere. 


3. Scattering of plane waves by a 
compressible sphere 


a) Plane progressive waves 

Let a sphere of radius a be placed freely at the 
origin of a co-ordinate system. The incident velo- 
city potential ©; of the plane progressive waves 
propagating in the direction of z-axis or 4 = 0 
is expressed by the following expansion 


PD; — eilot—kz) 


= elot ¥ (2n + 1) (— 


n=0 


i)” Jn(kv) Pr(cos8), (16) 
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where k is the wave number, J, the spherical 
Bessel function of the order n and P,, the Legendre 
function of the order n. 

The velocity potential outside the sphere, ©, 
may be written as 


D1 = 0, + ®, (17) 
where ®, is the velocity potential of the scattered 
waves and may be expressed as follows: 
oo (18) 
PD, = elt Y (2n + 1) (—i)” Ay HO(k v) Pp(cos 9) , 
n=0 


H" being the spherical Hankel function of the 
second kind, A, a constant to be determined by 
the boundary conditions. 
The velocity potential inside the sphere, 2*, 
may have the following expression: 
fe (19) 
O* — et) (2n + 1) (—i)” Ba In(k* v) Pn(cos 6), 
n=0 


where k* is the wave number inside the sphere — 
the mark * is used hereafter to indicate quanti- 
ties concerning the sphere—, By, a constant to 
be determined together with A, by the following 
boundary conditions, 


* 
se = = at y (20) 
and 0091 =0* @* at v=a, (21) 


oo and o* being respectively the density of the 
outer and inner medium. A, and B, are found 
as follows: . (22) 
re Ak Sn(k*a) J) (ka) —k* J) (k*a) In(k a) 

" k* J) (k*a) H2(ka) —AkIn(k*a) H®’ (ka) ’ 


ik 
(ka)? [k* Jp (k*a) HO) (ka) —Ak Jp(k*a) HY’(ka)] ” 
A= 0*|0o. 
If (ka)? <1, (k*a)2?<1 and A= O((ka)?}, 
Bo takes the following expression 
¥ 3i 
~ —(k*a)2ka+ i[32—(k*a)?] © 


where 


Bo (24) 


This is the case of a bubble suspended in liquid 
and Bo is proportional to the amplitude of the 
pulsation at the surface. Eq. (24) shows the 
resonance in the pulsation, which occurs at 


(k*a)2 = 3 A, 


the resonance frequency f, for an ideal gas bubble 
being given by 


(25) 


(26) 
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where y is the ratio of the specific heats of the 
gas, po the hydrostatic pressure. The calculations 
in the next section show, for plane progressive 
waves, a prominent peak in the radiation pressure 
corresponding to resonance in the pulsation of 


the bubble. 


b) Plane stationary waves 


If the centre of the sphere be at a distance h 
from a nodal plane, the incident velocity poten- 
tial referred to the equilibrium position of the 
centre as origin may be written 


DP; = ei l[etk +h) 4 e-ik E+h)] 


= eat (2 n+ 1) (—i)” On Jnlk v) P,,(cos 9), (27) 


where bn = (—1)melth + e-tkh, (28) 


The scattered velocity potential, ®,, and the 
velocity potential inside the sphere, ®*, may be 
written 


P, = clo! Ee n+1)(-i)” Cy H®(k ») Pp (cos), (29) 


D* — eit) (2n+ 1) (-i)” Dn In(k* v) Pn(cos6). (30) 
n=0 
Cy, and D, are determined by the boundary 
conditions given by (20) and (21) as follows: 
Cn os An On» Dy = By On. (31), (32) 


4. The radiation pressure on a 
compressible sphere 


It is convenient to express P,,, Pug, Puy, ue. 
Po, in terms of ©*, and by using the relations 


* 
= at vod, (33) 
09 BD; = o* o* NE ain (34) 
oD, ss, a" a 

00 ap = 0" aR AG = (35) 


together with the substitution cos) =p, we 
have from (12), (13), (14) and (15) the following 


expressions 


+1 
, [(aet\? 
«Pu, = naryo | ( Oy ) & Le du>, 


—l 


(36) 


hed 


o@m* 2 
(Pua = Catton {(F—) w—r)dw>, (87) 


—i 


+1 


(Papas) =<—2aaigo | ( 


—1 


op* 


ov 


‘op* 
a 
¥ a 


yaa \ Of , 


x (l— 2) du), (38) 
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oe Ns 


iY ee ena 


68 f (6*)?_, udu. (39) 


It is also convenient to write the solution | 


obtained in Section 3 in the form 


Re(D*),-. = (2n +1) MnPp(cos6), (40) 

Re ( ae ag ‘ =i (n+) K, Pteoshl.. (40) 
My = Re[e#!(—i)” By Jn(k*a)], (42) 

Kn = Re[et!(—i)n B, k* Ji,(k*a)]. (43) 
Substituting (40) or (41) into (36), (37), (38) or 


(39), we have the expression of the radiation 
pressure in terms of M,’s and Ky’s. 


<P> = CP.) + CP) 4 Pe) Cee eee 
UP. = <2 weg y2 (noel) (Ke Roa eae) 
i (46) 


GPig 20 AZ00 ys 2n(n+1)(n+ 2) <Mn Mnid, 
< Py, ue> =—2 24/1 09 L2 (n+ 1) (n+ 2)<Kn Mnii>+ 
+ 22aloo Y 2n(n+1)<Mn Kn+1>, (47) 
n=0 


(Po) =—2ark?a?2? 90 2 (n+ 1) Mn Myst. (48) 


If M, and Ky, are written in the form 


Mn = Min coswt + Men sinwt, (49) 


Kn = Ki,n coswmt + Ke,n sinwt, (50) 


we have 
2 <Mn Mnsi> = My, n Mi, n41 + Me2,n Me, n+1, (51) 
2 <Kn Mansi> = Kin Mi, n+ + Ko,n M2,n41 (52) 


and so on. 
Writing B, (in the plane progressive wave 
field) in the form 


if i 
Bn ka i Di, n—iDan z (33) 
where (54) 
D1,n=k*aJi(k*a) Jn(ka) —Aka Jn(k*a) Jak a), 
(55) 
Dz n=k*a Ji (k*a) Nn(k a) —AkaJ,(k*a) Np (ka) 


and Ny is the spherical Neumann function, and 
substituting (42), (43) with (53) into (44), (45), 
(46), (47) and (48), it will be proved that 
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QP = 200 Xx 

© (n+ 1)(Di,n De,n41—De2,nDi, n+}? 
x @0 2) D3 ,\|D? D3 

n=0 (Di nt talt i? nti t 2, n+] 


Hence, the radiation pressure on a compressible 
sphere in a plane progressive wave field is always 
positive, i.e., in the direction of wave propa- 
gation. 


. (56) 


5. Radiation pressure on small spheres, 
[(ka)2, (k*a)2<1], of which A= 0*/09 = O[1) 


This is the case of a small liquid sphere 
suspended in another liquid. If (ka)? and (k*a)? 
are negligibly small compared to unity M,’s 
and K,’s, for which n= 3, will contribute 
nothing to the final result. 

Consequently we have from (43), (44), (45) 
and (46) 

(P) = —2 700 [2 a2 (1+ 24) <Ko Ki> + 
+ 4a?2(A—1) (Ky Ke) + 2k? a2 72 (Mp My) + 
+124(1—A) <M Mo2)}. (57) 


a) Plane progressive waves 

Substituting M,’s and K,’s into (57), we ob- 
tain the result 

<P> = 227 0 (k a)§ F(A, 0) 

I 

2% 


(58) 


= 7 a2-4(ka)*— F(A, 0), (59) 
where o6 = c*/c=k/k*, I(=oock?/2) is the 


intensity of the incident waves and (60) 


1 142A\2 2 
(l+2a)2 (2 3102 +5 Oe | 
F(A,o) is the density-compressibility factor, 
which reduces to Kine’s density factor when 


o-oo. It is positive and the radiation pressure 
is in the direction of wave propagation. 


F(A,c) = 


b) Plane stationary waves 


Substituting M,’s and K,’s into (57), we 
obtain 


<P> = 427 00(k a)? sin2 kh F(A, 0) 
=za*4kaFsin2kh F(A, 0), 


(61) 
(62) 
where E(= 09k?) is the mean total energy- 


density in the stationary wave field and F(A,o) 
the density-compressibility factor given by 


4+ [2(A—1)/3] 1 
eae 3102” 


which approaches KiNnc’s density factor when 
o-oo. If F(A, a) is positive, the radiation 


F(A, 6) = (63) 
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wee 
—> 2-070, 
A+ [2(—D)/3] 1 


Fig. 1. The curve F(A, c) = Ea) 


27, 310 


and the results of experiments. 

(1) Toluene in water, (2) Benzene in water, (3) 
Nitrobenzene in water, (4) Mercury in water, (5) 
Ortho toluidine in water, (6) CCl, in 2.5 % Na,CO, 
aq. solution. 

N and L in the figure show that the accumulation 
was observed respectively at the nodes and at the 
loops. 


pressure urges the sphere away towards the loops 
and if negative towards the nodes. 

The curve F(A,o)=0 is shown in Fig. 1. On 
the right of the curve, F(A, c) is the positive and 
the spheres are to be forced away towards the 
loops and on the left towards the nodes. 

The points indicated by 1, 2, 3 and 4 are from 
the experiments by K. SOLLNER and C. Bonpy 
on emulsions, while the points 5 and 6 from our 
data (unpublished). According to Kinc’s theory, 
the critical density ratio is 0.40 and denser 
particles than 0.40 x 99 accumulate at the loops 
and vice versa. SOLLNER and Bonpy found an 
experimental rule that particles denser than the 
medium of dispersion accumulate at the loops 
and vice versa. The figure shows the best agree- 
ment between the experiments and the present 


calculations (Eq. (62) and (63)). 


6. Radiation pressure on bubbles 
[(ka)®, (k*a)?<1, A= Ol(ka)*}] 
a) Plane progressive waves 


Putting M,,’s and K,’s into (44), (45), (46), (47) 
and (48), we have 


ee > «Pu» K Pas ues <Po> (64) 
and 
KFS = Fue» (65) 
= —2 2a? 09 Ko,0 K2,1, (66) 
x 2meoar(k*a)& 
~ o2 (k*a)6 + [3 A— (k*a)2]2 ” ap 
4 x a? (k*a)4 I (68) 


~ e(o2 (k*a) + [3A —(k*a)?]?) * 


where I(=ock?/2) is the intensity of the inci- 
dent waves. 
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The frequency dependence of the radiation 
pressure given by Eq. (67) is shown in Fig. 2 for 
air bubbles in water (0°C, 1 atm). It has a 
prominent peak at the pulsation resonance, 


k*a=/3 A, and takes a minimum value at 


k*a~ 3. 
0.04 
0.03 + A 
ae 
Vv 
| 0.02 
0.01 a 
0 2 4 


— k*a/V32 


Fig. 2. The radiation pressure on air bubbles given by Eq. 
(67), at 0°C, 1 atm, the peak value being taken as 
unity. 


The ratio y of the magnitude of the radiation 
pressure on the resonant bubble calculated from 
Eq. (68) to that calculated from Kinc’s formula 
is found to be 

apni epyp se 1h 
BCP me 180 


Let f, and Af be respectively the resonance 
frequency and the width between two frequen- 
cies, at each of which the radiation pressure falls 
to a half of its peak value and the pulsation 
amplitude to 1/]/2 of its peak value. If the 
sharpness Qprp of the < P>-frequency curve is 
defined by the ratio f,/Af, we have (smallness of 
Af|f, being assumed) 


(69) 


pa 


cL 1 


D™ oV3a 


Ore = hr (70) 


r aS 
Af 
where Q, is “Q of the pulsation mode” and D 
the logarithmic decrement of the resonant 
bubble. 

In virtue of Eq. (70), we can determine the 
logarithmic decrement of the resonant bubbles 
from a study of the radiation pressure. 

For air bubbles in water at 0°C, 1 atm, we 
obtain from the known value of 2 and o 


y~8x10% Qy x70. (71) 


The experiments of the authors and Eq. (68) 
are found qualitatively in good agreement. 


Qn = 


172 


Quantitatively, however, the theory gives for 
both v and Q, values 7---10 times too great, if the 
effect of heat conduction is not taken into 
account. 


Calculations, in which the heat conduction’ 


effect in bubbles is taken into account, and 
results of an experiment will be given in another 


paper. 
b) Plane stationary waves 


Putting K’s into (65), we have 


<P> = —2. 2 0% 092 (Ko Ki), 


= —2207 09 Ki,0 Ki,1, (72) 
= —4 709 sin2 kh F(A, o, k*a), (73) 
or 
A Find kh F(A, 0, k* 14 
cere sin (A, o, a), ( ) 


where E is the mean total energy density in the 
plane stationary wave field and F(/,o,k*a) is 
given as 
o (k*a)[3 A— (k*a)?] 
F *q) = 

(A, 0, ka) = Sa pea)6 + [3A — (Rta) 


The frequency dependence of F(A, o,k* a) is 
shown in Fig. 3 for air bubbles in water (0°C, 
1 atm). If k*a> 34, F(4, 0, k*a) is negative 
and the radiation pressure is in the direction from 
the nodes towards the loops and if k*a <//3A, 
towards the nodes. At a given frequency, smal- 
ler bubbles than the resonant should be forced 
away towards the nodes and greater ones to- 
wards the loops. On the resonant bubbles, the 
radiation pressure in the plane stationary field 
should vanish. 

Observations by the authors, which will be 
given in another paper, and the above-stated 
result of the calculation are qualitatively in good 
agreement. 

0.6 


_ (75) 


F (4,0, ka) 


- 0.4 


~0.6 


Fig. 3. F(A,o,k*a)- k*a curve for air bubbles in water, 
at 0°C, 1 atm. 
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Appendix 
In the calculation of the radiation pressures 
on small spheres use was made of the following 
expressions for M,’s and K;,’s. 
a) Small spheres [(ka)2, (k*a)? <1], of which A =O\1). 


Plane progressive waves: 


(k*a)? (k*a)*ka(Ao2 —1) 
aKyjp=— cosmt + 922 
xX sinw t, 
(k.a)4(1—A) ka ; 
ak,= 31+ 2A2 cos@ t+ 7 og sing t, 
2 (ka)? 4 (ka)? (1—A) 
tha 30+3A) 18s ee 
X sino t, 
pik (k*a)2ka(Ao2—]) .. 
Mo=- cosot 3 72 sing t, 
ka)*(1—A) 
M, “3(1+ 202 cosmt + I 97 not 
Jr (k a)? 2(ka)*(1—A) . 
Mz=> 3(2 432) cosm@ t + 71352. (2 +302 sinw t. 
Plane stationary waves: 
* \2 
Pee eee a 2 coskh coswm t + 
(k* a)*k a (Ao? —1) 4 
972 2 coskhsinwt, 
ak sit Oi ink h cosw t + 
1= apy) sin S@ 
(ka)*(1—A) ,.. : 
os “31h 2d? (1+ 202 2sinkhsinwt, 
As 2 (ka)? 
a Ke = —3 9 gqy 1 One nceerae 


A (kay? (14) 
135(24 342. 


2 coskh sino t, 


Mico nee cue eee 


A 
(k*a)2k a (Ao? —1) : 
- 372 2coskhsinwt, 
ia 2h oe er t 
1= yy sin cos@mt + 
(k a)4 (1—A) . , 
30 (1+ 2A2 2A) 2 sinkh sinw t, 
(ka)? 
1) ee 
2 CRED ae 


2 (ka)? (1—) 


135 (2+ 3 22 y. coskhsinw Be 
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b) Small spheres [(ka)?, (k*a)? <1], of which 
A= O((k a)?}. 


Plane progressive waves: 


1 


aKo= Teah(era's Sa— (rape | (a) * 
x [3 A—(k*a)?] cosw t — (k*a)4 (ka) sing t}, 
4 
Dis Ma coswt+kasinot, 
kia)2 ? 
Arce = cos@ t + sing t, 
Mo d 


~ (ka)? (k*a)!+ [34 —(k*azp2 * 


x (3 [3 A—(k*a)?] coswt + 3 (k*a)2kasino t}, 


(k a)* : 
Mi =< coswt+ kasinwt, 


3 
ka)? (ka)? . 
Migr Cee 
2 rR coswm t+ 270 sing f. 
Plane stationary waves: 
* )2 
eee (k*a)? 2 coskh 


(k a)? (k*a)*+ [3 A— (k*a)?]? 


x (—[3 A — (k*a)?] cosw t — (k*a)? kasino t}, 
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LIS: 


aK, =—ka2sinkhcoswt + © (ka)tsink h sinw t, 
(ka)? 


aKe=———.— 2 coskhcoswt + 


(ka)? ‘ 
cae LP 2 coskhsino t, 


_ 2coskh 
(k a)? (k* a)4 + [3 A — (k* a)?]? 
x (3[3A—(k* a)?] coswt + 3(k*a)2kasino t}, 


Mo= 


x 


DAY 
M,=—ka2sinkhcoswt + CO" gsinkhsinot, 
2 
M,=—*) 2 coskh cosa t + 
(ka)? 
+ 270 2coskhsina t. 


(Received 12th January, 1955.) 
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ACOUSTIC RADIATION PRESSURE ON BUBBLES 
AND THEIR LOGARITHMIC DECREMENT 


by K.Yostoxa, Y. KAWASIMA and H. Hirano 


Institute of Scientific and Industrial Research, Osaka University 


Summary 


Radiation pressure on bubbles in a plane progressive sound wave field is calculated taking the 
effect of heat conduction into account. At 70 and 120 ke/s in water, measurements were made 
on H, and O, bubbles. The result of the experiment is in good agreement with theory, showing 
a prominent peak corresponding to the pulsation resonance. Logarithmic decrements of the 
resonant bubbles deduced from the data on radiation pressure are also in good agreement with 
the theory of Prrrem and Sanryosat. In the case of plane stationary waves, bubbles smaller than 
the resonant are found to accumulate at the nodes and larger ones at the loops, as expected from 
the authors’ calculations, while the radiation pressure on the resonant bubble vanishes. 


Sommaire 


On a calculé, en tenant compte de l’influence de la conduction thermique, la pression de radia- 
tion sur des bulles se trouvant dans un champ d’ondes sonores progressives planes. Les mesures 
ont porté sur des bulles de H, et de O, dans de l’eau, et pour les fréquences 70 et 120 kHz. Les 
résultats expérimentaux et théoriques sont en bon accord, et mettent en évidence un maximum 
accusé correspondant a la résonance pulsatoire. Les décréments logarithmiques des bulles en 
résonance, décréments déduits des résultats de mesure de la pression de radiation, sont également 
conformes a la théorie de PrrreM et SANEYOsHI. On a constaté, en bon accord avec les résultats 
des caleuls des auteurs, que, dans le cas d’ondes stationnaires planes, des bulles plus petites que 
les bulles en résonance s’accumulent aux nceuds, tandis que les bulles plus grandes vont aux 
ventres; par pilleurs, la pression de radiation est nulle sur les bulles en résonance. 
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Zusammenfassung 


RADIATION PRESSURE ON BUBBLES 
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Der auf Gasblasen in einer ebenen fortschreitenden Schallwelle wirkende Schallstrahlungsdruck 
wurde unter Beriicksichtigung der Warmeleitung berechnet. Messungen an H,- und O,-Blasen 
in Wasser bei 70 und 120 kHz, wobei man einen Spitzenwert des Schallstrahlungsdrucks bei der 
radialsymmetrischen Resonanz erhielt, ergaben gute Ubereinstimmung mit der Theorie. Auch 
die aus den Werten des Strahlungsdrucks berechneten logarithmischen Dekremente stimmten 
gut mit der Theorie von PrRteM und SANEyosuHI iiberein. In ebenen stehenden Wellen sammelten 
sich Blasen, die kleiner als die Resonanzblasen waren, an den Schwingungsknoten, die gréBeren 
in den Schwingungsbauchen, wie es von den Autoren auf Grund ihrer Berechnungen erwartet 
wurde, wihrend bei der Resonanz der Strahlungsdruck verschwand. 


1. Introduction 


Calculations of the radiation pressure acting 
on bubbles, which are freely suspended in a plane 
progressive sound wave field in a perfect fluid, 
are given in Section 2 taking the effect of heat 
conduction into account. 

In Section 3, the procedures and the results of 
measurements are given of the radiation pres- 
sures on Ho and Og bubbles in water. 

Of the logarithmic decrements of the resonant 
bubbles, the only data so far are by E. MEYER 
and K.Tamm™ [1]. They found rather larger 
values (especially in the ultrasonic range) than 
those indicated by the theories of H. Prriem [2] 
and J. SANEYOsHI [3], who calculated the decre- 
ment taking the effect of heat conduction into 
account. The values of the logarithmic decre- 
ment obtained from the data on radiation 
pressure, as well as the radii of the resonant 
bubbles, are in good agreement with the theories 
of PrrieM and SANEYOSHI. 


2. Radiation pressure on small bubbles, taking the 
effect of heat conduction into account 


In the present paper, the same notations are 
used as in the authors’ preceding paper [4] 
entitled “Acoustic radiation pressure on a com- 
pressible sphere”’. 

Let a bubble of radius a be placed freely in a 
perfect fluid, in which the train of the plane pro- 
gressive sound waves is propagating. Referring 
to the equilibrium position of the centre as 
origin, the incident velocity potential ®; and the 
scattered velocity potential ®, may be expressed 


D; a eilot—kz) — ery, (2n olin 1) (—i)” * 
n=0 (1) 


< In(k v) Pn (cos8), 


@, = elo)’ (2n + 1) (—i)” An H?) (kv) Pn(cos8). (2) 
n=0 

It is assumed that compared to the wave- 

length of the sound, the radius of the bubble is 

so small that (ka)?, (k*a)2<1 and that the 

motion of the bubble corresponding to n > 2 in 


the order of spherical harmonics contributes 
nothing to the final result. Even if the effect of 
heat conduction is introduced, the following 
relations may hold as in the non-heat-conducting 
bubbles: - 
Pax > eye CPay ues «Po 


<P> = «Pu 
= —2200-2a?-<Ko Ky), (3) 


and 


where K’s are given by 
—u, = Ko + 3 Kj cos8, (4) 
Ko = Ki,0 cosmt + Ko osinwt, (5) 
K, = Kj,1 coswmt+ Kz 1 sinwt, (6) 
u, being the radial velocity at the surface, and 
2 (Ko Ki> = Ky,0 Ki,1+ Ke,0 Ke,1- (7) 


The expression for K; may not be altered by 
the introduction of the heat conduction effect, 
since the vibration corresponding to n= 1 is 
independent of the pulsation of the sphere. 
Consequently we have 


Ky, 1=— k (ka)3/3, Ko, 1=> k. (8), (9) 


Taking the effect of heat conduction into 
account, PrrieM gave the following formula for 
the motion corresponding to n= 0 


greta tee 
» Ana? dt 
we iwa 3 (y—1)k* fe 
3 y p* 2 
x 2 (s* a coshs* a — sinhs* a) I"] dp*, 


iw 0* cha? 
(10) 


where V is the volume of the sphere, dp* the 
excess pressure at the surface, 


«* — thermal conductivity, 
cp = specific heat at constant pressure, 
ct = specific heat at constant volume, 
ae ai 
y = ole, 
* __ |/2%@ &p 
s* = /——., 
K 
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i (11) 
k(sa+1) 


coshs*a-«*(s*a—tanhs*a) + sinhs*a-«(sa+1)° 


(The quantities marked with and without * are 
for the gas and the surrounding medium respec- 
tively.) 

It is convenient to write the relation (10) in 
the f 

Dap se —w oo u, = (G+ iF) dp*. (12) 
Since «*(s*acoshs*a)<x(sa-+ 1)sinhs*a and 
tanhs*a ~1 in the present experiments, F' and 
G take the following forms, 


(13) 
(k*a)? 3 y—1 sinh2/*a—sin2l*a 
F= 5 ; 
: 3A alle l*a ee | 
au Mee 3 (y—1) 

tle aan Gore 

( sinh21*a + sin2l*a i 14 

cosh2 1*a — cos2l*a ae) aid) 
where 

A= 0* /e0, 
a = % esaye 2 oh 


(eee 


Determining Ag in (2) so as to satisfy the 
relation (12), we obtain 


1 a F(a F—1) + (aG)? 


0 ae kaa)? 4 (¢F—1)2’ (15) 
a aG+ka(aF)? 
-Sadgeiraicri+@Fr—i: “° 
which show the resonance to occur when 
aF=1. (17) 


Substituting (8), (9), (15) and (16) into (3), the 


radiation pressure is found as 


(ay? + =2 8° 
<P) = Pro we *) ; (18) 
nel +(1-Sp) 
| 4 ak 
or 
ka-aG 
CET 


I being the intensity of the incident sound. 
The formulae show that, corresponding to the 
pulsation resonance, the radiation pressure takes 
the maximum value 


(20) 
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where (ka), is given by Eq. (17) as 
(ka), =o (k*a), 


eRe ; an lenis 
=oV8i (1+ 5-4, 1 sinh2/*a—sin21 


(21) 


a ye 
a, 


If, as in the present experiment, the last term 
in the bracket in (13) or (21) is negligible, we 


have 


a cosh21*a — cos2I* 


N D, 
<P) mae = aa Ge S MAL Bo a st Bd (22) 
7 
<P) maz = “5—> (23) 


<P max being the limiting value of ¢(P)maz, 
when the thermal conductivity «* +0. Di, and 
D, are respectively the logarithmic decrement 
of the resonant bubble arising from the heat 
conduction and from the scattering of the sound 
and by the sum D,+ Dip is given the logarith- 
mic decrement of the resonant bubble due to the 
combined effect. D;, and D, have respectively 
the forms 


Din =2G/F, 
D,=cV32. 


(24) 
(25) 


If A(ka) is the width between two values of 
ka, at each of which ¢ P) falls to half of the peak 
value, and the sharpness Qpp of the <P)-(ka) 
curve is defined by the ratio (k a),/A(k a), we ob- 
tain from (18) or (19) 


1 Alka) | D+ Dy 1 
Qrp (ka), I Or 


where Qpis the “Q of the pulsation of the bubble”. 
In virtue of the last relation, we can determine 
the logarithmic decrement of the resonant bubble 
from the experiment on the radiation pressure. 


(26) 


3. Procedures and results of the experiment 


a) Apparatus and procedures 


The experimental arrangement is schemati- 
cally shown in Fig. 1. A water tank of 50 x 35 x 
35 em? is lined with sheets of Acrilite, on which 
a number of rectangular cones of the same 
material, having a bottom area of 1 cm? and a 
height of 1.6---2.4 em, are closely affixed. T is a 
piezoelectric composite transducer, serving as the 
sound source. 

Bubbles of various dimensions are formed by 
electrolysis on a small platinum electrode P. 
When bubbles come up to the axial region of the 
sound beam, the switch in the transmitter circuit 
is closed and the train of sound waves is sent out. 
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Fig. 1. The experimental arrangement. 


The path of the bubble is then deflected and 
when the sound intensity is not too great, the 
motion is found to reach a stationary state in 
a short time. 

Photographs of the tracks of the bubbles were 
taken by intermittent illumination, an example 
being shown in Fig. 2. When the sound intensity 
is too great, the tracks show zigzag figures. 


=“ 
eed 


1 
t 
% 
*) 
' 
14 
ee | 
i 
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Fig. 2. Tracks of H, bubbles in water, the time interval 
of illumination being of 0.28 s. The radius of the 
bubble on the left: 0.052 mm. 70 ke/s approx. 
50 uW/cem*. White spots are microscopic flow indi- 
cators. 


Within the Fraunhofer region, the observa- 
tions were made where the sound field as cal- 
culated was sufficiently close to plane progressive. 

The upward velocity vy; in the absence of the 
sound and the horizontal component v2 of the 
terminal velocity in the sound field are measured 
on photographic prints. 

From the upward velocity y;, the radius of the 
bubble was determined assuming Stokes’ law for 
the motion of spheres in viscous fluid. According 
to the work of ALLEN and others, departures 
from the Stokes velocity are negligibly small if 
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the radius is less than 0.06 mm in the tempera- 
ture range of the present experiments in water. 

From the horizontal component velocity ve, 
the radiation pressure was determined by the 
formula CP enaae 
where 7 is the viscosity coefficient of water. 

In order to distinguish the motion of the 
bubble relative to the fluid from that due to the 
direct flow caused in water by the second-order 
effect of the sound, a small amount of microscopic 
precipitate was suspended in water as flow indi- 
cator. The observation of the Tyndall light 
showed throughout the present experiments no 
trace of motion of these precipitates, even when 
the bubbles were in vigorous motion. 


b) Results obtained in the plane progressive wave 
field 


An example of the result obtained at 120 ke/s 
on Oz bubbles in water is shown in Fig. 3. The 
radius of the resonant bubble is found to be 
0.027 mm and the sharpness Qrp 9.0, from which 
we obtain the logarithmic decrement of the- 
resonant bubble as 0.35. 


0.15 
mm/s 
= 
® 
S 0.10 + 
3 
' 0.05 
0 0.02 0,03 004mm 
——> Radius 


Fig. 3. The radiation pressure on O, bubbles in water, at 
120 ke/s, 9°C..The electric input to the transducer 
390 uW/em?. 


The values obtained for the radius of the 
resonant bubbles are listed in the last column of 
Table I. The radii a, in the table are calculated 
by the formula k*a=|/32. The radii ag are 
obtained from the equation aF = 1 and multi- 
plied by /1+ 2T/(ap), T and p being respee- 
tively the surface tension of water and the 
prevailing hydrostatic pressure. The discrepan- 
cies observed with Hz bubbles may be attributed 
to their sensitiveness 
water. 


to contamination from 
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Table I 
The radius of the resonant bubble 


Radius of the resonant bubble 
in mm 


Frequency 
in ke/s 


calc. dy 


experimental 


Table II 
The logarithmic decrement of the resonant bubble 


The log. dec. of the resonant 


Frequency Cas bubble 
in ke/s 
calculated experimental 
120 
70 
70 


The values obtained for the logarithmic decre- 
ment of the resonant bubble are listed in the last 
column of Table II. The values of the decrement 
in the third column are calculated from known 
data by Dj, + D, using the expressions (24) and 
(25). 

The agreement between the calculated and the 
observed values is considered as satisfactory. The 
latter are also indicated in Fig. 4. 


0,8 


—— Logarithmic decrement 


0 50 70 100 
——> Frequency 


120 kc/s 


Fig. 4. The logarithmic decrement of resonant bubbles. 
Solid curve: calculated for O,. Broken curve: 
calculated for Hy. 

e obtained by Meyer and Tamm for Og. 
© obtained by the present experiment for O,. 
x obtained by the present experiment for H,. 


c) Results obtained in the plane stationary wave 


field 


A plane reflector, made of sponge-gum sheet, 
was placed parallel to the wave front in the beam 
of the sound. O» bubbles of various dimensions 
were formed by electrolysis on a platinum wire 
electrode at the bottom of the tank. When 
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bubbles came up to the region in front of the 
reflector, the sound waves (of 69 ke/s) were 
transmitted. As is shown in Fig.5 tracks of 
bubbles were photographed by intermittent 
illumination. 

The reflection at the reflector surface was 
examined using a sonic interferometer and was 
found to be practically the same as that at the 
free surface. The reflector surface thus forms an 
antinodal plane, from which successive nodal and 
antinodal planes are located on photographic 
prints. 


~<i— loop 
<_— 

~<— loop 
~<i— node 
= __ loop 
<= node 


t sponge-gum reflector 


Fig. 5. The tracks of O, bubbles in the plane stationary 
field in water. Smaller bubbles than the resonant 
are accumulating at the nodes and those larger at 
the loops. The track a shows a penetrating resonant 
bubble. The bubble d is produced by the coalescence 
of b and ce. 


Examination of the tracks of bubbles showed 
clearly three groups of them (Fig. 5). The radii 
and the motional characters of some selected 
bubbles are listed in Table III. The bubbles of 
radii distinctly smaller than the radius of 
resonance accumulate at the nodes, those larger 
at the loops, quite in accordance with the authors’ 
calculations. 

The positions of the nodes and the loops are 
likely to have no effect upon the motion of the 
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Table III 


The motional characters of O, bubbles in the plane station- 
ary field, 69 ke/s, in water. The radius of the resonant 
bubble: 0.046 mm 


Remark 


Remark 


to node to node 
0.028 | to node penetrating to loop 
0.032 | to node penetrating to loop 
0.038 | to node to loop 


resonant bubbles, which is characterized by the 
penetration of the latter through the nodes and 
the loops away from the source, the velocities 
observed being much slower than in the pro- 
gressive field of equal mean energy density. This 
fact may be explained as follows. 

If the sound field is perfectly plane stationary, 
the radiation pressure on the resonant bubble 
may vanish as in the theory. If any difference 
exists, as is usual, between the intensities of the 
component waves, the resonant bubble will 
proceed through the nodes and the loops in the 
direction of wave propagation of the stronger 
component. This view is supported by the fact 
that the penetrating motion is more vigorous in 
front of the somewhat weak reflector (metallic 
plate) than in front of the nearly perfect reflector 
(sponge-gum sheet). 

In ultrasonic fields of low frequencies and of 
low intensities, k*a of bubbles, produced by 
cavitation in the vicinity of the nodes, may be of 
small value and the chances of coagulation few. 
In such cases, the accumulation of the bubbles 
may be observed at the nodes only, as reported 
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by R.W. Boye, G.B.Taytor and D.K. 
FRoman [5]. 

In ultrasonic fields of higher frequencies and 
of higher intensities, the lives of bubbles smaller 
than the resonant may be extremely short be- 
cause of the frequent chances of coagulation. In 
such cases, the accumulation of the bubbles may 
be observed mainly at the loops, as reported by 
O. A. ANGERER and F. Srrersi [6] as well as 
G. E. GoLpMAN and G. R. Rico [7]. 

It may further be expected that abrupt 
transitions of the bubbles may take place from 
the nodes towards the loops when the radii are 
increasing over that of.resonance. In the authors’ 
experiment at 460 ke/s and at the intensity of 
0.1---0.2 W/cm2, the accumulations at the loops 
were distinctly observable in a test-tube, while it 
was very difficult, but not impossible, to discover 
the faint and unstable nodal accumulations and 
the transitions of bubbles towards the loops from 
the nodes. 

A part of the present work was aided by the 
Scientific Research Expenditure Fund from the 


Department of Education. 
(Received 12 January, 1955.) 
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Buchbesprechung 


O. Wimes, Lirmbekimpfung, Grundlagen und 
Ubersicht. Gildeverlag, Alfeld/Leine 1954, 176 Seiten, 
40 Abb., DIN A5, kartoniert DM 6,—. 


Das vorliegende Biichlein stellt den ersten Band einer vom 
Deutschen Arbeitsring fiir Larmbekampfung herausgegebe- 
nen Schriftenreihe dar. Es enthialt eine Auswahl von Vor- 
tragen, welche von namhaften Fachleuten auf dem Ham- 
burger Anti-Liirm-KongreB im Jahre 1953 gehalten wur- 
den. 

Zur Kennzeichnung des Inhaltes seien kurz die behandel- 
ten Themen mitgeteilt: EinfluB des Larms auf den mensch- 
lichen Organismus; Aufgaben und Aussichten der Betriebs- 
larmbekimpfung; MeBtechnische Grundlagen der Lirm- 
bekampfung; Die Liastigkeit des Lirms; Schallschutz in 
Wohnbauten; Schallschutz auf Motorschiffen; Gerdusch- 
starke und -messungen an StraBenfahrzeugen; Lirm und 
Verkehrssicherheit; Ursachen des Lirms im Eisenbahn- 


verkehr und Erfolge in der Lirmabwehr; Lirmbekimpfung 
in der Luftfahrt; Rechtsgrundlagen der Larmbekampfung; 
Administrative MaBnahmen gegen den Larm. 

Es liegt hier ein Ubersichtswerk vor, das sich in knapper 
Form mit den medizinischen, technischen und juristischen 
Fragen der Auswirkung und Bekimpfung des Larms be- 
faBt. So bietet es vor allem dem interessierten Laien, der 
MaB8nahmen zur Larmminderung ergreifen will, einen guten 
Uberblick iiber den Stand der Entwicklung auf diesem 
Gebiet. Aber auch der Fachmann wird diesen zusammen- 
fassenden Bericht begriiBen, insbesondere wird den mit den 
technischen Problemen bereits vertrauten Leser die Be- 
handlung rechtlicher und medizinisch-psychologischer Pro- 
bleme interessieren. Durch die Mitteilung von Diskussions- 
beitragen wird die Lebendigkeit der Darstellung noch er- 
hoéht, und es werden jeweils die aktuellsten Fragen beriihrt. 


Wolfgang Westphal 
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Summary 

The effect of atmospheric turbulence on the overall sound intensity received downwind of a 
source at different distances from it and at various heights above the ground was investigated. 
The results revealed the occurrence of a zone of maximum intensity followed by another of mini- 
mum intensity at a higher level. The effect of “‘eddies” formed in the turbulent atmosphere on the 
diffraction pattern of the source was studied. 


Sommaire 

On a étudié influence d’une turbulence de l’atmosphére sur l’intensité sonore totale recue 
en aval d’une source, a différentes distances de celle-ci, et 4 différentes hauteurs au-dessus du 
sol. Il existe, d’aprés les résultats obtenus, une zone d’intensité maxima, suivie d’une autre zone 
d’intensité minima, 4 un niveau plus élevé. On a étudié l’influence des tourbillons de l’atmosphére 
turbulente sur le diagramme de diffraction de la source. 


Zusammenfassung 

Die Wirkung atmosphiarischer Turbulenzen auf die Schallausbreitung in Windrichtung wurde 
untersucht, indem man die Schallintensitat in verschiedenen Abstanden von der Schallquelle und 
in verschiedenen Héhen vom Erdboden bestimmte. Dabei ergab sich eine Zone maximaler Inten- 
sitat, der in gréBerer Hohe eine Zone minimaler Intensitat folgte. Weiterhin wurde der EinfluB 
von ,,Wirbeln“‘, die sich in der turbulenten Atmosphire bilden, auf die Richtcharakteristik der 


Schallquelle untersucht. 


1. Introduction 


A great deal of experimental research work in 
recent years has been devoted to the study of 
fluctuations in sound intensity received from a 
source of constant output in the atmosphere. It 
is clear that the additional attenuation which an 
open-air sound suffers as compared to an indoor 
one is due to turbulence in the atmosphere 
resulting in local fluctuations in temperature, 
pressure and wind speed. Results of observations 
are usually stated in one of two ways; either as 
mean values of sound intensity recorded during 
intervals extending over a few minutes (during 
which the average meteorological conditions are 
noted) or as instantaneous values of both sound 
and meteorological conditions. Each type of 
observations has its merits. The former method 
is used when a study of overall intensity is needed 
for giving an average picture to be used for 
sound ranging or signalling over different terrains 
at different times of the day. Instantaneous record- 
ing is used to give a picture of correlated events 
occurring simultaneously. The early investiga- 
tors on the propagation of sound in the atmos- 
phere adopted the large-scale recording method 
and explained their results in terms of the mean 
values of the prevailing meteorological conditions. 


BEHRENS [1] gave a comprehensive summary 
of the theory of sound diffraction caused by wind 
and temperature gradients. KOMERALL [2] dealt 
with the effect of temperature gradients on the 
sound path. Eacieson [3] and others studied 
the effect of temperature, wind and humidity 
on sound over distances ranging up to 200 me- 
ters. 

Working with both types of observations, 
RicHarpson [4] has recently studied the fine 
structure of a turbulent atmosphere in relation 
to the propagation of sound therein. He investi- 
gated also the overall sound attenuation at differ- 
ent levels from the ground for two types of 
terrain of different natures. The Authors [5] in a 
previous investigation studied the effect of wind, 
lapse rate and scale of turbulence on sound 
propagation over a desert terrain. 


The work described herewith comprises two 
more investigations and is complementary to 
the previous work. The next section deals with 
the estimation of overall sound attenuation at 
various distances downwind and at various 
heights above a sandy terrain. The concluding 
section deals with the effect of “eddy diffu- 
sion” in a turbulent atmosphere on sound pro- 
pagation. 
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2. Overall attenuation of sound 


The overall loss in the intensity of sound in 
the atmosphere is mainly due to three factors, 


namely, (a) the spreading which takes place in ' 


accordance with the inverse square law of 
distance r, amounting to 20 logr dB, (b) the 
molecular absorption which is controlled by the 
physical and chemical constitution of the air and 
(c) the scattering caused by the reflection, refrac- 
tion and diffraction effects and enhanced by 
gustiness of wind, eddy formation, adiabatic 
lapse rate and roughness of terrain. Factors (a) 
and (b) are now fairly well established. For the 
study of the third factor the authors have adopted 
the large-scale averaging method, in view of the 
fact that the nature of such a problem dealing 
with turbulence necessitates that averages only 
be subjected to analysis. Moreover, the recording 
of the meteorological conditions had to be done 
by instruments which yield mean values over 
finite time intervals. 

The experiments were conducted on summer 
days at a site North-East of Cairo over a desert 
terrain which was almost free from obstructions. 
The source of sound used was a loudspeaker 
actuated by a beat frequency oscillator (General 
Electric type 913) producing a pure tone. The 

6 2 
| { 
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output of the oscillator was led into an amplifier 
and then to a matched loudspeaker unit. The 
sound was received by a small crystal microphone 
whose output was passed to a sound level meter 
(General Radio type 759) having a slow response, 
so that a time average of intensity could be 
recorded. To reduce the background noise, the 
signal received was passed through a sharply 
tuned wave analyser acting as a filter (General 
Radio type 736A). The intensity close to the 
source was about 95 dB. The source was placed 
at different heights from the ground ranging 
from | to 6 feet. The microphone was placed at 
the same heightas the source and the distance 
between the two varied from 20 to 100 feet. 
Frequencies between 1 and 5 ke/s were used. 
The mean value of the wind speed was recorded 
for each observation at the same height as the 
sound beam. The difference of temperature was 
also recorded between heights of 1 and 4 feet 
from the ground. 

Results for overall sound intensity are exhibited 
in Fig. 1. The general behaviour of most of the 
curves differs from that in RicHaRDson’s results 
(cf. Fig. 6 of [6]) which had been carried out 
under different meteorological conditions, namely 
under less gusty weather, low humidity and 
small lapse rate. Few of our graphs are similar 

to his, possibly be- 
© cause of temporary 
local uniformity at 
: the time of measure- 
requency 

1000cs. = ment. 
v Summer desert 
conditions are char- 


ok 
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acterised by high 
lapse rate especially 
in the hours round 
noon. This results in 
the formation of 
strong eddy convec- 
tion currents which 


t 2000 c/s 
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give this characteris- 
tic behaviour to the 
curves. Such beha- 
viour can be explain- 
ed on the basis that 
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Fig. 1. The overall sound 
intensity for the 
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large sound dissipation takes place at ground level 
due to absorption by sound while in the layers 
just above, the intensity increases rather quickly 
and is further assisted by the sound reflected from 
the sand which is almost in phase with the direct 
beam.At a higher level, the ground dissipation dis- 
appears, but phase difference between the direct 
and the reflected beams starts to develop result- 
ing in a reduction of the intensity. Moreover, 
deviation of path takes place due to refraction 
in these layers; the actual path being approxi- 
mately parabolic. Wind gradients tend to 
decrease the overall sound intensity picked-up, 
shown by the concave part of the curve and 
these effects are more marked the greater the 
lapse rate. In the figures shown, the minima 
occurred at a height of 2 to 4 feet above ground 
level, when the lapse rate was at its maximum. 
At a higher level still, the phase difference 
changes sign once more and intensification is 
produced; at the same time, the upward refrac- 
tion decreases due to the gradual fall in lapse 
rate with height. If observations are continued 
at higher levels still, another minimum, but less 
marked, is likely to appear. 


3. The effects of eddy diffusion on sound 
propagation 

Turbulence in general is a direct consequence 
of high lapse rate, and implies a state in which 
the instantaneous wind velocities exhibit irre- 
gular and random fluctuations. Such fluctuations 
take the form of “eddies” having a size ranging 
from a few inches up to several feet. 

In order to demonstrate the role of eddy 
diffusion on the propagation of sound, measure- 
ments of intensity were taken in horizontal 
planes across the sound beam downwind, using 
a small condenser microphone and frequencies 
ranging from 1 to 5 ke/s. The common height of 
source and receiver was varied from 2 to 6 feet 
above ground, while the distance of the line of 
cross traverse of the microphone from the source 
was maintained at 5 feet. Sample results for a 
frequency of 3000 c/s are plotted in Fig. 2 in 
which the theoretical diffraction pattern corres- 
ponding to a baffled piston source is drawn as a 
thick broken line for comparison. 

It is clearly seen from the figure that the out- 
door observations, while approximately main- 
taining the general diffraction pattern, exhibit 
certain departures from it, such as: 

(a) The maximum intensities of the main central 
beam fall below that recorded under exactly 
similar conditions in an anechoic chamber 
free from turbulence. 
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Fig. 2. Effect of “eddy” diffusion on sound diffraction in 
turbulent atmosphere; 

2 ft above ground, 

4 {t above ground, 

6 ft above ground, 

assumed theoretical distribution. 


(b) The maximum intensities of the diffracted 
beam are definitely higher than their corres- 
ponding values in the theoretical curve. 
The minimum intensities in the experimental 
curves are not zero as required by theory. 
The bearing of the outdoor minima from line 
ahead is greater than that calculated from 
the theoretical relation sinf = 0.61 A4/R 
where / is the wave-length and R is the 
radius of the source, and 6 is the angle to the 
first minimum. 

At high levels above the ground the “wings” 
tend to broaden out. 

Summing up all these effects, we may conclude 
that atmospheric turbulence results in general 
loss in overall intensity of the direct beam and 
smudging of peaks of the diffracted “wings”, thus 
indicating that the phenomenon is essentially 
one of scattering rather than one of simple 
refraction. 

The authors hope to investigate the effect of 
negative lapse rate on sound propagation and to 
extend their observations to cover higher levels 
than hitherto explored. 
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MESSUNGEN DER REFLEXION VON BIEGEWELLEN 
AN QUERSCHNITTSSPRUNGEN AUF STABEN 


von Muciono* 


III. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 

Die vorliegende Arbeit dient der Untersuchung der Ausbreitung von Biegewellen in Bau- 
elementen. In der eindimensionalen Vereinfachung, namlich fiir Biegewellen auf Stangen, wurden 
die an einfachen und doppelten Querschnittsspriingen auftretenden Reflexionen theoretisch 
berechnet und experimentell gemessen. Experiment und Theorie zeigen befriedigende Uberein- 
stimmung. 


Summary 

The aim of this work was the investigation of the propagation of flexural waves in building 
structures. In the simplified one-dimensional case, i. e. for flexural waves in rods, the reflection 
coefficient for single and double changes in cross-section has been calculated and measured. 
Experiment and theory show good agreement. 


Sommaire 

Le présent travail est destiné 4 faciliter l’étude de la propagation des ondes de flexion dans les 
éléments de construction. On a calculé théoriquement et mesuré expérimentalement, dans le cas 
unidimensionnel simple, c’est-a-dire pour des ondes de flexion dans des verges, les réflexions qui 
ont lieu a des discontinuités de section simple ou double. Les résultats expérimentaux et théoriques 


sont en accord satisfaisant. 


1. Einleitung 


Die Fortleitung von Biegewellen auf Wanden 
und anderen Bauelementen ist fiir die Bau- 
akustik von grofem Interesse, da die Schall- 
dammung zwischen zwei Raumen stets, unab- 
hangig von der Starke der Trennwand, durch die 
Schallabstrahlung der zu Biegeschwingungen er- 
regten gemeinsamen Seitenwainde eine obere 
Grenze findet. Dariiber hinaus spielen auch die 
auf der Trennwand selbst angeregten Biege- 
wellen (Koinzidenzeffekt) eine wichtige Rolle. 
Fiir die Planung von Bauten ist es deshalb wich- 
tig, Klarheit zu schaffen iiber die Ausbreitungs- 
bedingungen fiir Biegewellen in Bauelementen 
allgemein und iiber die Reflexion von Biege- 
wellen an Isoliermaterialien, Ecken und Quer- 
schnittsspriingen im besonderen. Theoretisch 
wurden die hierbei auftretenden Probleme von 
L. CREMER [1]  eingehend behandelt. Experi- 
mentelle Untersuchungen machten EXNER und 
BoumeE [2], die an Staben die Schalldammung 
von zwischengeschalteten Isoliermaterialien ma- 
Ben und Kurrze, Tamm und VoceEt [3], die die 
Ausbreitungsverhaltnisse an Ecken mit und ohne 
zwischengeschaltete Isolierstoffe wiederum an 
Hand von Stabmodellen studierten. 


* Djokjakarta, Indonesien. 


Die vorliegende Arbeit beschaftigt sich mit der 
durch einfache und doppelte Querschnitts- 
spriinge hervorgerufenen Schalldammung bzw. 
Reflexion fiir Biegewellen. Untersucht wird wie 
in den vorerwahnten Arbeiten der eindimensionale 
Fall, indem die Messungen an Staben ausgefihrt 
werden. Gemessen wurde an folgenden Anord- 
nungen: 


a) Einfacher Querschnittssprung (Abb. 1) 


In dem skizzierten Stab, dessen Querschnitt an 
der Stelle x= 0 bei gleichbleibender Breite von 
dem Wert h, seiner Hohe auf den Wert hz diskon- 


tinuierlich  springt, 

werden am linken ‘ | f 

Ende Biegewellen(in = (1)_— jy ha (2) 
der Zeichenebene) i | } 
erregt, die sich tiber 
den Querschnitts- 0 


sprung hinweg nach Abb. 1. 

rechts fortpflanzen. 

Das rechte Stabende ist (siehe unter ,,MeBver- 
fahren‘) reflexionsfrei abgeschlossen. Die an 
der Sprungstelle reflektierten Biegewellen werden 


gemessen. 
b) Doppelter Querschnittssprung (Abb. 2) 


Hier tritt der gleiche Querschnittssprung 
«x =he/h; = he/hg doppelt auf, und zwar wurde 
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sowohl der Fall einer Querschnittserweiterung 
als auch der einer entsprechenden Querschnitts- 
verengung untersucht, wobei sich die Frage- 
stellung auf die Abhangigkeit des Reflexions- 
faktors von « und dem Abstand / zwischen den 
beiden Querschnittsspriingen bezieht. 


Abb. 2. 


2. Theorie 
Fall a) 


Die Wellengleichung fiir Biegeschwingungen 
lautet!: (1) 


Be 
Dy, Wz, Mz, F,), 


22 r 

ape Ou» Wz, Mz, Fy) = — a on Taman (6 
wobei x die Fortpflanzungsrichtung, v, die 
Schnelle in y-Richtung, wz die Winkelgeschwin- 
digkeit um die z-Richtung, M, das Biegemoment 
um die z-Achse, F', die Kraft in y-Richtung, B die 
Biegesteifigkeit und m die Masse pro Langen- 
einheit in x-Richtung darstellen. 

Fiir sinusférmige Erregung kénnen wir setzen 


cee OO Te) 


Setzt man dies in die pr Sap Chung ein, 


so erhalt man 
BkiA=mo2 


und wegen w/k = cp (k= 22/Ap, cz = Fortpflan- 
zungsgeschwindigkeit fiir Biegewellen) 
cp =)/B/m Vo: (2) 


Im vorliegenden Fall gilt mit dem Elastizitats- 
modul E, der Dichte 9 und der Stabhéhe h: 


4 
BO weet 
ca= Vis Vho, (3) 
AOS Ae 
Setzen wir als Lésung fiir die Differential- 


gleichung @,(x,t) = i,(x)e)°', so erhalten wir 


die gewéhnliche Differentialgleichung: 
B d By(x) _ 


m dx 


«? by (x) (5) 


1 Der Querstrich iiber den Zeichen bedeutet, daB es sich 
um komplexe GréSen handelt. 
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mit der Lésung (6) 


by(x) = Aye“ik2 + Apetika + Bye-kx +4 Bo etkz, 


Bei dem in Abb. 1 dargestellten diskontinuierli- 
chen Ubergang miissen an der Ubergangsstelle 
die folgenden vier Beziehungen gelten: 

Wz1 = Wz2, 


Ma a Mz. 


Dy1 — Dy2 > 


2 re (7) 
By = Fy, 


Bei Biegewellenerregung am linken Stabende 
wird an der Ubergangsstelle ein Teil der Biege- 
wellen reflektiert und der Rest durchgelassen: 


Dyl = e—ihax+ R etikc R’et+ kya, (3 
dye = De-ih2 4 D'e—he, 

Substitution der Gleichungen (8) in (7) liefert, 
unter Beriicksichtigung der Beziehungen 6, = 


dy/dt, #, = diy/dx, M, = —(B/jw) diw,/dx und 


Fy = —dM,/dx, die vier Gleichungen: 
heh pp 41, 
jR+R4+jKD+ KD'=j, 
—R+R+yD—yD'=1, 0) 


aR Rak yD ak Doi, 
worin K und y die folgenden Bedeutungen haben: 
= ke/ky = ms By /m, By = A\/22, 
y =k} By/ki By =\/mz By/m, B, . 


Fur den Fall, da®B an der Ubergangsstelle nur 
eine Durchmesser- aber keine Materialanderung 
auftritt, ist 


K=o-¥2 und yp=a?. 

Als Lésung der Gleichungen (9) ergibt sich 

“4 2 yp (1 K?) jK(1 y)? (10) 
K(it yt + 2y(1+K% 


Ps) 


Der durchgelassene Bruchteil der Energie ist 
d=1— |Ri2. Definieren wir als Reduktionsfaktor 


[R] = 10 log(1/d), so gilt (11) 


5/4 4. y-3/4 4 43/4 a ~2 
[RI =10 os AU a al a 


po + x2 + 14 oll? + - + 2 


und haben damit die Lésung fiir den ersten Fall 
in der Form 


Fall 6b) 
Betrachten wir zunachst beliebige Querschnitts- 


spriinge (Abb. 3), so gilt: 
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Abb. 3. 


Vyl == eikiaz + Ry etikiax + Cy ethia, 
Dy2 ot As edhe at 4 Ro etikee + Co eke ® 4. ro ek. (al), 
(12) 


Zweimalige Anwendung der Gleichungen (12) fiir 
x = 0 und fiir x =/ liefert die acht Gleichungen: 


14. R,+ €, =Ao+ Ret Co+ Copel, 
Ag e-iksl + Ro etilal + Eye! + C5 = Az t+ Co, 
ky (—j+jRi+ G)= ka (—j Ae + jRo—C, + 
+ Cy e—*l), 
ko (—j Ap eS! + j Rp et tel — Cy ee + C)) 
= k3(—j As — Cs), 


Vy3 — Az e—iks(%—1) + C3 e—ks(x—l) , 


B, ky Bd oe bo 
one g a FG) (13) 
= i. 2 (—Ay — Ro + Co 4+ Cy e-*l) 
@ 
as a (— Ap elk! — Ro et ial 4 Cy el + C5) 
k3 B 
a POSS (ag ae Ge 
jo 


kB; duet oo ee ee ee 
x "re Cite jite— Ge + Cae), 
kp B, | Ao eniksl — 4 Ry etikzt — Co e-kal 4 C! 
Bias (+ j Ag e Ske! — j Ro et ike! — Cy e- ks! + Cp) 
ke Bae ee 
— jo (j Az — C3). 


Setzen wixijetzt h;—h,, also) K =k, /k; = hy /k, 
=o? und y= ki B,/k? By) Bike Boo 
und vernachlassigen wir alle Terme mit e~*:!, 
so ergeben sich schlieBlich die acht nachsten 
Gleichungen mit den acht Unbekannten R,, 
C,, Az, Ro, C2, Cs, Az, C3, aus denen Rj, die 
Amplitude der reflektierten Welle, ausgerechnet 
werden kann. Die genannte Vernachlassigung ist 
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Abb. 4. Blockschaltbild der MeBanordnung. 
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Abb. 5. Das benutzte Kérperschallmikrophon. 


erlaubt, da bei den Versuchen / so gewahlt wurde, 
daB 3a<kol ist und damit e—*:! <10-4. 
Exner und BoumeE [2] zeigen auBerdem, dab das 
Nahfeld bereits in einem Abstand von 4 von der 
Ubergangsstelle nicht mehr mefbar ist. Im vor- 
liegenden Fall wurde 1 = 1,52 gewahlt: 
R, an GG. deere =—l, 
jR,+C,4+ K €,+jK A,—jK R25j, 
—R,+ C,—ypCo,+y4o+yR.=1, 
—jRi+ C+ KyCo— jKy42.+jKyR.=~-j, 
(14) 
+ eIhl4, + elk! Ro + Ci —A3—C3=0, 
—j Ke-ik:l 4, + j Keikz! Ro + KC. ot j As aia C3=0, 
—peik14,—pelks! Ro + p C3 + A43—C3=0, 
+jK pe thl 4, — 
j K peikl Ro+ Ky C,—j A3+ C3=0. 


Als Lésung ergibt sich fiir Ri: 


ee 


j N2 e+ ik! — (L + j M)(—M + jL) eit! 
2 N(M cosk21— L sinkg 1) 


(2 L M coskz1— (L2— M2 + N2) sinkol] 


j ((L?— M2 — N2) coskgl + 2L Msinkg1] (3) 
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mit L=2y(1—K2), 
M=— K(1—y)?, 
N=K(1+y)?+2y(1+ K%, 
a =he/hy=he/hs, po2, K=o-1/?, 


Aus der Bedingung: M coskz:!1— L sinkzl=0 
ergibt sich die Wellenlange, fiir die der Reflexions- 
faktor verschwindet. Man sieht, daB diese Wellen- 
lange nicht nur von der Lange | des Ubergangs- 
stiickes abhangt, sondérn auch vom Querschnitts- 


verhaltnis «; Ry = 0 fir kol = arctg M/L. 


3. MeSverfahren 


Die benutzte MeBanordnung ist im Block- 
schaltbild in Abb. 4 wiedergegeben. Die benutz- 
ten Stabe sind 3---4 m lange, relativ verlustarme 
Duraluminiumstangen, die an diinnen Faden auf- 
gehangt wurden. Die Faden wurden gegen die 
Stative und diese wiederum gegen den FuBboden 
mit Weichgummi isoliert. Fir die Messungen 
wurde ein kleines Kristallmikrophon (Beschleu- 
nigungsmesser nach H. OBERsST) verwendet, das 
durch seine geringe Masse die Wellenausbreitung 
auf der Stange nicht beeintrachtigt (Abb. 5). 
Das Gehause des Mikrophons besteht aus Elek- 
tron; der Kristall ist ein einseitig eingespannter 
ADP-Bieger, der feuchtigkeitsunempfindlich ist. 

Empfangsseitig wurde zur Erhéhung der Mef- 
empfindlichkeit ein Oktav-BandpaB und ein RC- 
Resonanzverstarker benutzt. Die Frequenzmes- 
sung erfolgte mit Hilfe von Lissajous-Figuren 
unter Benutzung einer Quarzuhr-Normalfre- 
quenz von 500 Hz auf dem Bildschirm eines 
Oszillographen. 

Fiir die Messungen muf das rechte, vom Sen- 
der abgewandte Stangenende unendlich lang, 
d. h. reflexionsfrei abgeschlossen sein. Ein solcher 
reflexionsfreier AbschluB kann z. B. dadurch er- 
zielt werden, daf man das Stangenende keil- 
férmig schneidet und in einen mit Sand [2] oder 
Sagemehl gefiillten Trog einfiihrt; man kann 
auch das keilférmige Stangenende mit einem Ent- 
dréhnungsmittel (wir haben ein Gemisch von 
Sand und Gummilésung im Gewichtsverhaltnis 
4:1 genommen) umgeben [4]. Beide Abschliisse 
erwiesen sich als ausreichend reflexionsfrei, d. h., 
die Reflexionsfaktoren waren bei den benutzten 
Frequenzen kleiner als 10%. 

Als Sender wurde ein kraftiges permanent- 
dynamisches System benutzt, dessen Schwing- 
spule fest mit der Stange verbunden war. Die 
Schwingungsrichtung der angeregten Biegewellen 
war stets die y-Richtung, d. h. die Richtung, in 
der die Stabmafe verandert wurden (Abb. 4). 
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Gemessen wurde immer der Reflexionsfaktor 
des Querschnittssprunges, und zwar durch Ab- 
tasten des senderseitigen Stangenteiles mit dem 
Beschleunigungsempfanger. Aus den Minima und 
Maxima der hier auftretenden stehenden Wellen 
ergibt sich der Reflexionsfaktor R nach der be- 
kannten Formel 

il a. 
Pe Se 


mit d’ = Vmin/Vmaz, dem Verhaltnis der Schnellen 
bzw. Beschleunigungen in Minimum und Maxi- 
mum. 


4. Me8ergebnisse 
Fall a) 


Mit einer Aluminiumstange entsprechend Abb. 
1, bei der der rechte Stangenteil (hg) quadrati- 
schen Querschnitt hat, wurden Messungen fir 
verschiedene Héhen h, des linken, senderseitigen 
Stangenteiles gemacht. hg betrug 40 mm und hy 
wurde variiert entsprechend « = hg/hy = 2; 5 
und 11,4. Fir die einzelnen Falle wurde der 
Reduktionsfaktor [R] = 10 log(1/d) berechnet, 
wobei d= 1— R2 den Bruchteil der einfallenden 
Energie darstellt, der durchgelassen wird. Es 
ergaben sich die folgenden Mefwerte: 


ie 5 oo LA 
Hz [R] Hz -fR] Hz ({[R] 
500 = 0,26 500 4,1 500 ~=-8,0 
1000 = 0,22 1000 ~— 3,3 1000 8,5 
1500 = 0,48 1500 = 3,5 
2000 =0,46 2000R 3a): 


Bei « = 11,4 wurde nur bis 1000 Hz gemessen, 
da bei héheren Frequenzen wegen der kleinen 
Schallgeschwindigkeit auf der diinnen Stange die 
Wellenlange nicht mehr groB gegen die Mabe des 
Kristallmikrophons ist. 


— Reduktionsfaktor [RJ 


1 2 5 10 20 
——~ ha /hy 
Abb. 6. Reduktionsfaktor [R] als Funktion des Quer- 
schnittsverhaltnisses x. Die ausgezogene Kurve ist 
berechnet, die Punkte geben die MeBresultate 
wieder. 
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Die theoretische Berechnung von [R] liefert: 


eae 
2 ae 
Bien tee 
114. 483 


Der Reduktionsfaktor [R] ist also theoretisch 
unabhangig von der Wellenlange. Die bei den 
Messungen beobachteten Abweichungen sind aus 
Abb. 6 ersichtlich, in der die ausgezogene Kurve 
die theoretischen Werte und die Punkte die 
Resultate der Messungen wiedergeben. 


Fall b) 
Die Wellenlange der Biegewellen ist gegeben 


durch i is 
E h 
ip 2a) ie ‘ 


In diesem Fall? haben wir eine sehr verlustfreie 
Messingstange benutzt, bei der sich experimen- 
tell ergab: 


Ap = 198 Vhjo. 


Wahlt man nun hy = h3 = 10mm und hg = 
15 mm (Abb. 2), so mu die Lange | des Zwi- 
schenstiickes so gewahlt werden, dab kol23a 
ist. Nimmt man | = 750 mm, dann ist 


bei 360 Hz: kel = 3,00a, 
bei 821 Hz: kel = 4,402. 


Die fiir diesen Fall gemessenen Werte sind zu- 
sammen mit der theoretischen Kurve in Abb. 7a 
wiedergegeben. 

In einem zweiten Versuch wurde eine Quer- 
schnittsverengung mit « = 0,52 gemessen. Mit 
1 = 750 mm ist 


bei 406 Hz: 
bei 690 Hz: 


kol = 5,30 a, 
kol = 6,90 a. 


2 Die ersten Versuche in Gottingen ergaben eine zwar 
kleine, aber bestimmt wesentliche Diskrepanz zwischen 
Theorie und Experiment, nimlich eine Verschiebung der 
experimentellen gegeniiber den theoretischen Kurven in 
der Frequenzrichtung. Der Grund dafiir diirfte sein, daB 
die verwendete Stange aus drei zusammengenieteten Teilen 
bestand. Bei nichtidealer mechanischer Verbindung sind 
Abweichungen zu erwarten, die sich als eine Frequenzver- 
lagerung auswirken. Wir haben daher die Experimente in 
Delft (Prof. KostEN) mit einer Stange wiederholt, deren 
Teile durch Hartlétung miteinander verbunden waren. Die 
beschriebenen Resultate beziehen sich auf die Messungen 
in Delft; sie weisen die erwahnte Diskrepanz tatsachlich 
nicht mehr auf. 
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Die Resultate zeigt Abb. 7b zusammen mit den 
berechneten Werten. 

Es kann daher als SchluBfolgerung gesagt 
werden, daB die MeBergebnisse mit der Theorie 
recht gut ibereinstimmen. 


+ aati — Mites Seranbamrorsmcaes (al 
| 


052) ae 


06 T 


607" 


a at kal 


64n 


687 


Abb. 7. Reflexionsfaktor R, eines doppelten Querschnitts- 
sprunges in Abhingigkeit von der Frequenz 
(a) Querschnittserweiterung, 
(b) Querschnittsverengung; 

berechnete Kurve, °°° Mefpunkte. 


Herrn Prof. Dr. E. Meyer, Géttingen, bin ich 
fiir den zur Verfiigung gestellten Arbeitsplatz in 
seinem Institut, fiir die Themenstellung und 
wertvolle Ratschlage zu Dank verpflichtet. 
AuBerdem danke ich Herrn Prof. Dr. C. W. 
Kosten, Delft, fiir seine Anregung zu der Ar- 
beitsrichtung. Herrn Dr. G. Kurrze habe ich fiir 
seine Unterstiitzung bei der Abfassung der 
Arbeit zu danken. 

Die Arbeit wurde teilweise mit Mitteln des 
Bundeswohnungsbauministeriums durchgefihrt. 

(Eingegangen am 10. September 1954.) 
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ATTENUATION OF THE (1,0) AND (2,0) MODES IN 
RECTANGULAR DUCTS 
by A. M. GHABRIAL* 
King’s College, Newcastle upon Tyne, England 


Summary 


The attenuation constant of the (1,0) and the (2,0) modes of a rectangular duct have been 
studied in the frequency range 1180---3400 c/s. A rectangular duct, confining a constant volume 
between a source and a reflector, is constrained to move between fixed rails. The various modes 
are excited by the ‘hole’ method and the sound field within the duct is surveyed by a probe 
microphone small enough not to disturb the sound field. The standing-wave ratio is measured by 
a resonance peak method which has the advantage of avoiding the need to delineate the complete 
standing-wave pattern along the tube and it is suitable for determining small attenuations which 
is the case in this work. The experimental attenuation constants for the (1,0) and the (2,0) modes 
agree with the theory of BrEatry within 1 %. 


Sommaire 

On a déterminé dans la gamme de fréquences 1180---3400 Hz la constante d’affaiblissement 
des modes d’un conduit de section rectangulaire. Ce dernier se déplace entre des rails fixes tout 
en conservant entre source et réflecteur un volume constant. L’excitation utilise la méthode du 
«trou», et on explore le champ au moyen d’une sonde microphonique. On mesure le taux d’ondes 
stationnaires par la méthode du maximum de résonance qui permet de ne pas déterminer la 
totalité du diagramme le long du conduit et convient pour mesurer les petits affaiblissements. 
Les constantes d’affaiblissement trouvées expérimentalement pour les modes (1,0) et (2,0) con- 
cordent a 1 % prés avec celles déduites de la théorie de Beatty. 


Zusammenfassung 

Die Dampfungskonstante der Schwingungstypen (1,0) und (2,0) in einem rechteckigen Rohr 
wurde im Frequenzbereich 1180---3400 Hz gemessen. Um die hier auftretenden geringen Damp- 
fungen noch messen zu kénnen und die Abtastung des ganzen Schallfeldes zu vermeiden, wurde 
zur Bestimmung der Welligkeit nur die Resonanzspitze der Rohrresonanz mit einem Sonden- 
mikrophon ausgemessen. Die so ermittelten Dimpfungskonstanten stimmen innerhalb von 1 % 


mit der Theorie von Brarry iiberein. 


1. Introduction 


The attenuation of higher modes in tubes has 
been investigated theoretically by P. M. Morse 
[1] considering a small wall admittance. RocGERs 
[2] derived theoretically an approximate expres- 
sion for the attenuation constant in rectangular 
tubes having absorbing walls. Harric and 
LAMBERT [3] measured the attenuation of the 
(1,0) mode in a rigid rectangular tube having a 
slot along its length in the frequency range 
2400---4000 c/s and put forward an empirical 
formula for the attenuation in the gas due to 
viscosity and thermal conductivity. LAMBeErr [4| 
later derived an approximate theoretical expres- 
sion for this attenuation constant. Beatry [5| 
and SHaw [6] used CREMER’s [7] boundary layer 
method to get expressions similar to that of 
LAMBERT for the attenuation of higher modes in 
ducts. Suaw tested his theory for the (1,0) mode 
in rectangular ducts and in different gases. 


* Now at the University of Cairo. 


This brief survey indicates that little work has 
been done on the measurement of the attenuation 
of the higher modes in ducts, and even the 
measurements of LAMBERT and SHAW were for 
the (1,0) mode only. 

In this work the attenuation of the (2,0) mode 
as well as that of (1,0) mode is measured in air 
by an improved method. 


2. Apparatus 
2.1. The wave guide 


The wave guide is a rectangular tube 300 em 
long 17.2 cm wide, and 8.9 cm in height, made of 
hardwood 2.5 em thick. It is bolted by rectangu- 
lar wooden belts 5 em wide and 2.5 cm thick 
fixed 30 cm apart to keep the cross-sectional 
dimensions constant within 0.2%. By the aid 
of wooden guides, the tube is constrained to move 
along its axis by coupling it to a large slow- 
motion screw reading to 0.025 mm. These guides 
as well as the wooden binders a foot apart 
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Wooden 
frame 


- dial 


oo rs: 


Moving 
carriage 


Fig. 1. (a) Sliding carriage, (b) Vertical traverse for micro- 
phone. 


prevent flexural vibrations of the walls on the 
tube. 

A slot 0.5 em wide is cut across the middle of 
the tube (Fig. la). A wooden sliding plate 1.5 cm 
thick, having a spline which fits the slot is used 
to carry the microphone across the tube. Its 
length is so chosen that it always covers the 
slot completely whatever may be the position of 
the microphone M. In order to reduce any possi- 
ble vibration of the walls of the tube near the 
slot, this part of the tube is strengthened by a 
2 cm thick wooden board. 

A vertical moving device (Fig. lb) is fixed to 
the sliding carriage and carries the microphone M. 


2.2 The source 


The construction of the source is shown in 
Fig. 2. CFED is a brass uniform circular tube, 
1.5 cm diameter. In the limb DE is fitted an air- 
tight piston P, pierced by a 0.3 cm diameter hole. 
This piston acts as a source and is connected to 
an ear-phone L, used as a loudspeaker, by a hard 
rubber tubing R. The second limb CE is fitted 
with a moveable solid piston P», acting as a per- 
fect reflector. These two pistons can be moved 
by screws. 

The part EF of the tube is flush to the end 
plate AA of the main wave-guide, and two holes 
for the (1,0) mode or three holes for the (2,0) 
mode half a wave-length apart are symmetri- 
cally drilled through AA and EF. Although in 


theory it is advisable to make the spacing equal 


A. M. GHABRIAL: ATTENUATION IN RECTANGULAR DUCTS 


ACUSTICA 
Vol. 5 (1955) 


to a half wave-length, yet in practice it is found 
that this spacing is not very critical. A change of 
the resonance frequency by about 10% does not 
affect the condition. The positions of the holes 
are chosen to be at pressure nodes with respect 
to the resonating tube CDEF, this choice being 
most convenient as the impedance at a pressure 
node is nearly infinite. Thus the impedance of the 
end plate AA of the wave-guide can be consid- 
ered infinite irrespective of the presence of the 


holes. 


A ZZ F 


Reflector 
piston 


Fig. 2. Source end of guide. 


2.3. The detector 


A measurement of the pressure within the tube 
is carried out by a probe microphone. The probe 
is an injection needle of about 0.05 cm in dia- 
meter bent at a right angle at its far end to 
minimise its reflection effect on the sound field 
in the plane where the measurements are taken. 
It is connected to an ear-phone, operating as a 
microphone, by a small piece of rubber tubing 
to reduce the effect of the mechanical coupling 
between the tube and the microphone. There is 
no need to calibrate the probe microphone since 
it is sufficient to obtain relative measurements at 
a specified frequency. 


3. Sources of error 


3.1. Probe reflections 


Sound reflections from the probe arise either 
from the body of the tube or the probe itself, or 
from its input admittance. The part arising from 
the body of the tube is minimized by choosing a 
probe of the smallest possible dimensions with 
respect to the wave-length, and by bending the 
end of the probe so that the measurements can 
be taken on a plane other than that confining the 
body of the probe. 

Any obstacle which is small compared with the 
wave-length in the wave-guide behaves like a 
shunt admittance on the equivalent transmission 
line. This analogue enables us to derive the 
following simple expression for the coefficient of 
reflection of a probe: 
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_ 10 ~ (40 — Np) 
» no + (no + np) 2 


where 79 is the admittance at a point in the tube 
before the probe of input admittance 7» is 
introduced. 

If 7» is very small in comparison with 7, then 
the reflection from the probe is negligible, for (1) 
“pa Tp = —Np|2 No- (2) 
This condition can be fulfilled by choosing a non- 
resonant probe or else by confining the measure- 
ments around a pressure antinode where the 
admittance 79 is obviously large. 

_ In the present work these two conditions are 
fulfilled by choosing the length of the probe for 
any frequency to be non-resonant, and by 
making measurements near pressure antinodes. 

The following test was made to ascertain the 
effect of the probe on the sound field. A constant 
length sound interferometer used for the test is 
shown in Fig. 3. A and B are two circular tubes 
of internal diameters 6.0 and 5.8 cm respectively, 
sliding telescopically one inside the other; 5 is the 
source and R is a heavy and smooth reflector. 


Fig. 3. Test of sound reflections from microphone probe. 


In making the test the tube is first adjusted 
to resonance. The position and value of a certain 
minimum is recorded by a probe P, the probe M 
under investigation is inserted into the tube B 
at a distance of a quarter wave-length and the 
position and value of the previous minimum are 
again checked by the probe microphone P. No 
change could be detected. However if the probe 
M is inserted at a pressure node and the same 
procedure is repeated, a very small change of 
about one per cent in the value of the pressure 
minimum is observed. 


3.2. Vibration of the tube walls 


To test how far the attenuation is affected by 
wall vibration, the attenuation coefficient was 
measured in a preliminary tube without any 
‘binders’. The result was very high compared 
with the theoretical value and the discrepancy is 
attributed to wall vibrations, flexural and other- 
wise. To eliminate these the tube was gripped by 
thick wooden belts at distances of a foot apart 
and was constrained between two rigidly fixed 
rails. 
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3.3. Constancy of the cross-sectional dimensions 


The cut-off frequency of the wave guide for 
any higher mode depends on its cross-sectional 
dimensions, and variations of these dimensions 
will result in unequal spacing of the pressure 
antinodes in the guide. 

This source of error was reduced by the use of 
wooden belts to keep the cross-sectional dimen- 
sions of the tube constant to within about 0.5 % 
along the whole length. Moreover measurements 
have been confined to a part of the tube where the 
observed wave-length was constant to about 0.2%. 


3.4, Non-linearity of the acoustical system 


In the present work the relative response of the 
microphones P and M was compared with the 
output of the oscillator to eliminate the effect of 
changes in the electrical circuit. However it is 
necessary that the acoustic system should be 
linear and to check this the voltage fed into the 
loudspeaker was increased in steps of 1.5 volts 
and the response of the microphone measured. 
Perfect linearity existed up to 30 volts and the 
working voltage was always taken below 3 volts. 


4. Methods of measurement 


4.1. Procedure and theory 


The general procedure is to adjust the sound 
source to excite the required mode at the required 
frequency and then the tube to resonance as 
indicated by the maximum response of the micro- 
phone M at the reflector end. The presence of the 
pure mode is then checked by surveying the 
field across the tube. The wavelength within the 
duct is measured by locating the successive 
minima along it. 

Now before describing the method of meas- 
uring the attenuation constant due to the gas in 
the wave guide, it may be well to express the 
attenuation constant in terms of measurable 
quantities. 

According to P. M. Morse, for a solid termi- 
nation in the standing wave pattern of the plane 
wave, the magnitude of the excess pressure P, is 


| Px| = |Po| //cosh2 « x + cos2rx (4) 
where « is the attenuation constant in neper 
em! and r is the reflection coefficient at the 
termination from which is measured. 

If x2 <r? as in the present case, the standing 
wave ratio is 

oS Pinin / Paz = tanha x (5) 

tanh"! S. (6) 


or C= -~ 
x 
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If tanh! is plotted against x a straight line 
should be obtained whose slope gives «. 

The correctness of the suppositions made in 
deriving the equation appears justified from the 
linearity of the experimental relation between 
tanh1S and x (Fig.4 and 5). Moreover, the 
attenuation constant in the present work changes 
from about 140-10 up to 200-10 nepers 
em-!, while the wavelength within the guide 
changes from about 10 to 60 cm. 


4.2. Measurement of the standing-wave ratio 


The standing-wave ratio S can be measured by 
recording the sound pressure along the wave- 
guide with particular attention to the values at 
the nodes and antinodes, but in carrying out 
these observations, it is found that the values of 
the pressure minima are too small to measure. 

The alternative method is basically analogous 
to the curve-width method of determining the 
Q-factor of a tuned circuit. The width W of the 
trough in the standing wave pattern is measured. 
This width is defined as the distance between two 
points of equal pressures on both sides of a 
pressure antinode, at which the value of the 
pressure is K times that at the antinode. K may 
be chosen arbitrarily and according to conven- 
ience. In this way the standing-wave ratio S can 
be expressed in terms of K and W/A, where Ag is 
the guide wavelength. 

The pressure P, at a point near an antinode in 
a standing wave pattern is given by 


Pz, = P, (1+ 7) (7) 


where r, is the reflection coefficient at the point, 
and P, is the pressure in the wave progressing in 
the x-direction. Now 


ry = |r| 2 ev (8) 


where y=2a/Ag is the propagation constant. 
Here x will be the distance of the point in ques- 
tion from the antinode. 

In this treatment it is assumed that the ampli- 
tude of the reflection coefficient remains constant 
around the antinode, since the attenuation « as 
well as the distance x are small quantities. 

The magnitude of P, is given by 


| Pa|? =| P|? (1+ |r?|—2|r| cos2 yx). (9) 


If x —W/2 then by the definition of K it follows 


that 
; xe Lt Ir? —2Ir| cosy 


1+ |r|? —2 |r| 


(10) 


Now the standing-wave ratio is given by S = 
(1—|r|)/(1+ |r|). So relation (10) reduces to 


K2=1+ (S2—1) sin?y W/2 (11) 
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whence 

s siny W/2 sin (2 W'/Ag) 

~ VK2—cos2y W/2 K2—cos2y W/2- 
Me Be 


This equation enables the standing-wave ratio to 
be calculated from readily measurable qualities. 

To determine K, the sensitivity of the ampli- 
fier and the power of the oscillator are increased 
so that the meter indicates a considerable deflec- 
tion, about 30 scale divisions (readable with a 
certainty of about one quarter of a scale division) 
when the microphone is at an antinode. The 
other observations around the antinode are 
conducted in the usual way. In practice, it is 
preferable in most cases to plot a graph for the 
relative pressure against the position of the 
probe-microphone M and to calculate W and K 
therefrom. 


(12) 


5. Results and discussion 


Results have been taken at six different fre- 
quencies for each mode. For the (1,0) mode the 
frequencies range from 1180 up to 1200 c/s. For 
the (2,0) mode the frequencies are embraced in 
the interval 2220 to 2400 c/s. 

Figs. 4 and 5 represent the standing-wave 
ratio tanh 1S against the distance x from the 
reflector. The variation of the attenuation in the 


0 20 40 60 
—+ Distance from reflector 


80 cm 


Fig. 4. tanhS vy. reflector distance; (1,0) mode. 


0 20 40 
——~ Distance from reflector 


Fig. 5. tanh"'S y. reflector distance; (2,0) mode. 
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Fig. 6. Attenuation v. wave-length of plane wave; (1,0) 
mode. 
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Fig. 7. Attenuation v. wave-length of plane wave; (2,0) 
mode. 


duct with the plane wave length for the (1,0) and 
the (2,0) modes is shown in Figs. 6 and 7. These 
figures show the large increase of the attenuation 
of the higher modes in the duct as the frequency 
approaches the cut-off frequency of a particular 
mode. 

Following SHaw the expression for « can be 
written in the form 


a (4.7 + Ag +) Vf (13) 
: g 
where 
x Vab (2, 1 
Aya FEE + (yy AP (+ =) 
and 
es Wage 2 (15) 
c bee 
Here 
c¢ = sound velocity in the an nf 
A = wavelength oi eee See 
C, = specific heat of the gas at constant volume 


= ratio of the principal specific heats of the gas 
kinematic viscosity of the gas, 

= k/yCyo =thermal ‘diffusivity’ of the gas, 
where k is the coefficient of thermal con- 
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ductivity, and 9 is the mean density of 
the gas, 

z; = width of the guide, 

y1 = depth of the guide, 

f =} frequency. 


On insertion of accepted values of constants for 


air in eqs. (14) and (15); 


Ay = 2.375 - 10-6 em! s!/?, 
Ag = 2.249 - 10-6 cm! s!/?, 


By substituting these values in (13) the ‘theo- 
retical’ attenuation constant is obtained. 

Tables I and II show the agreement between 
the ‘theoretical’ attenuation constant calculated 
from the above expression (13) and the measured 
attenuation constant as found from Figs.4 and 5. 


Table I 
Values of attenuation for (1.0) mode 
(cut-off frequency = 967 c/s) 


Frequency theoretical experimental 
c/s nepers cm! nepers cm-* 
1180 140 - 10-5 143 - 10-5 
1240 134 - 10-5 135 - 10-5 
1320 130 - 10-° 130 - 10-5 
1395 LT stOK* 127 - 10-5 
1550 125 - 10-5 124 - 10-5 
1700 126 - 10-5 125 - 10-5 


Table II 
Values of attenuation for (2,0) mode 
(cut-off frequency = 1935 c/s) 


theoretical experimental 


Frequency 


c/s nepers cm-} nepers em-! 


6. Conclusion 


By using a narrow non-resonant probe tube 
with a bent end situated at a pressure antinode 
the disturbing effects of the probe in the acoustic 
field have been considerably minimised. This is 
particularly important in attenuation measure- 
ments. The technique employed to measure and 
calculate the standing-wave ratio from the width 
of a trough is especially suitable for small 
attenuations, as besides reducing reflection errors 
at the pressure maximum it avoids the uncer- 
tainty in determining, by the ordinary way, the 
magnitude of the pressure minimum. 

The attenuation constant measured for the 
(1,0) mode is in agreement with the theoretical 
expressions of Lampert, Beatry and SHAw. 
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The attenuation for the (2,0) mode is in con- 
formity with the theory to about 1% which is 
within the limits of experimental error. 
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EIN EXPERIMENT ZUR DEMONSTRATION VON 
KAVITATIONSBLASEN UND KAVITATIONSKEIMEN 
von W. Gira 
III. Physikalisches Institut der Universitat Gottingen 


Zusammenfassung 


Es wird ein einfaches Gerit beschrieben, mit dem sich die Umrisse von Kavitationsblasen in 
etwa 100facher VergréBerung projizieren und dadurch einem groBen Horerkreise sichtbar machen 
lassen. Die Blasen entstehen an einer Stahlmembran, die durch einen Druckstof angeregt wird. 
Kurze Zeit nach der StoSanregung wird die Kavitationsblase yon einem Hochspannungsfunken 
beleuchtet und dadurch projiziert. Bevor sie sich ausbildet, entstehen unter dem EinfluB des 
durch den Sto8® hervorgerufenen Unterdruckes viele kleine Blasen an im Wasser befindlichen 
Kavitationskeimen. 


Summary 


A simple apparatus for projecting the outlines of cavitation bubbles 100-fold for demonstration 
to a large audience is described. The bubbles are formed on a steel diaphragm excited by a pressure 
pulse. Shortly after the excitation, the cavitation bubbles are illuminated by a high-voltage spark 
and so projected. Before they grow there appear, under the influence of the negative pressure 
set up by the pulse, many small bubbles on the cavitation nuclei existing in the water. 


Sommaire 


On décrit un dispositif simple qui permet de projeter, avec un grossissement d’environ 100 fois, 
la silhouette des bulles de cavitation, de fagon 4 montrer ces bulles 4 un grand nombre de specta- 
teurs. Les bulles sont produites par une membrane en acier qui est excitée par une impulsion 
de pression. Un temps court aprés cette excitation, les bulles de cavitation sont éclairées par une 
étincelle haute tension et sont projetées sur un écran. Avant l’apparition d’une bulle de cavitation 
de nombreuses petites bulles se forment 4 partir des germes de cavitation existant dans l’eau, 


sous l’effet de la dépression créée par Vimpulsion. 


1. Einleitung 


Kavitation ist der Ubersetzung des Wortes ge- 
maf die Bildung von Hohlraumen. Und zwar 
sind hier stets Hohlraume in Fliissigkeiten, ins- 
besondere in Wasser, gemeint. Man kann sie er- 
zeugen, indem man das Wasser, oder kleine 
Bereiche darin, einer Zugbeanspruchung solcher 
Starke aussetzt, da die Zugkrafte die Krafte 
der molekularen Bindung iiberschreiten. Dann 
zerreiBt das Wasser, es bildet sich ein Hohlraum, 
eine Kavitationsblase. Auf diese Weise entsteht 
Kavitation in Diisen, an Schiffschrauben, an 
Wasserschallsendern usw. 

Allerdings brauchen fast niemals die zur Er- 
zeugung von Kavitationsblasen notwendigen 


Zugkrafte die Krafte der molekularen Bindung 
zu erreichen. Das hat seinen Grund darin, dai 
stets kleine Fremdkérper — etwa Staub — vor- 
handen sind, zu denen die Adhasionskraft ge- 
ringer ist, als die Kohasionskraft des Wassers. 
Hinzu kommt, daB solche Staubteilchen meistens 
in eine diinne Lufthaut eingehiillt sind. Natiir- 
lich setzt an derartigen Fremdkérpern die Kavi- 
tation leichter und eher ein als im freien Wasser. 
Sie stellen also Kavitationskeime dar. Als Keime 
sind auch Rauhigkeiten — haufig wieder mit 
Lufthauten und -einschliissen — an festen Ober- 
flachen anzusehen. 

Fir ein Demonstrationsexperiment sind die 
oben genannten Verfahren der Kavitations- 
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erzeugung deshalb nicht sonderlich geeignet, weil 
sie einen apparativen Aufwand erfordern, dem 
das Inventar einer fiir Vorlesungszwecke be- 
stimmten Sammlung physikalischer Gerate in 
der Regel nicht gewachsen ist. Vor allem aber ist 
ja eine wesentliche Forderung, die man an einen 
Demonstrationsversuch stellen sollte, die nach 
Einfachheit und Ubersichtlichkeit der verwende- 
ten Versuchsanordnung. Mit dem Ziel, dieser 
-Forderung méglichst zu entsprechen, wurde der 
im folgenden beschriebene Versuch zusammen- 
gestellt und verwendet, der es ermdglicht, die 
Umrisse von Kavitationsblasen durch Projektion 
einem grofen Hérerkreise sichtbar zu machen. 


2. Der Versuch 


In einer mit Wasser gefillten Kiivette ist eine 
Stahlmembran fest eingespannt. Zwischen ihr 
und einer in geringem Abstand dariiber befestig- 
ten metallischen Spitze (Abb. 1) wird ein Hoch- 
spannungsfunke ausgelést. Der dadurch im Was- 
ser itiber der Membran entstehende sehr heftige 
DruckstoB beult die Membran ein wenig aus und 
treibt das darunter befindliche Wasser nach 


plats : 
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so kann man, je nach dem zeitlichen Abstand 
zwischen StoB und Belichtung, die Kavitations- 
blase in einer der eben beschriebenen Phasen 
ihrer Lebensgeschichte sehen, bzw. projizieren. 


3. Die Versuchsanordnung 


Kine Skizze des elektrischen Teils zeigt Abb. 2. 
Es ist im wesentlichen die Macusche Verzége- 
rungsschaltung [5]. Eine Hochspannungsquelle H 


Ane , 
gia eel 


Abb. 2. Macusche Funkenverzégerungsschaltung. 


ladt die Kondensatoren Cy, Co, C3 auf, bis bei 
etwa 8000 V die Durchschlagspannung der hinter- 
einandergeschalteten Funkenstrecken F; und F» 
erreicht ist. F, ist eine Luftfunkenstrecke, Fs die 
unter Wasser befindliche StoBfunkenstrecke zwi- 
schen der Spitze S und der Membran M. Durch 
das Umschwingen des aus C,, Cz und L bestehen- 


@ © © 


Abb. 1. Skizze der StoBanordnung und des Kavitationsablaufs. S Metallstab, M Membran, F, Stob- 
funkenstrecke. Die Pfeile deuten die Bewegung des Wassers an. 


unten in der Richtung der Pfeile (Abb. la) fort. 
Wegen der der Membran eigenen Riickstellkraft 
wird deren Auslenkungsgeschwindigkeit kurz 
danach geringer. Das Wasser aber ist dank seiner 
Tragheit bestrebt, die ihm durch den Stof iiber- 
tragene Anfangsgeschwindigkeit beizubehalten, 
und wird sich, wenn der Stof heftig genug war, 
von der Membran losreiBen, sobald die Adhasions- 
kraft iiberschritten wird. Die Membran schleu- 
dert also das Wasser regelrecht von sich fort; es 
entsteht eine Kavitationsblase (Abb. 1b). 

Wiahrend sich im weiteren Verlauf die Grenz- 
flache von der Membran entfernt, strémt von 
den Seiten her das dort befindliche Wasser in den 
Hohlraum nach, solange, bis die Blase sich véllig 
abgeschniirt hat (Abb. lc,d,e). SchlieBlich implo- 
diert dann die Blase, wie es von anderen Beob- 
achtungen iiber Kavitationsvorginge her be- 
kannt ist [1]---[4]. 

Wird einen Augenblick nach der Stofaus- 
lésung das Gebiet unter der Membran mit einem 
sehr kurzzeitig brennenden Lichtblitz beleuchtet, 


den elektrischen Kreises kehren sich die Ladun- 
gen an den Platten von Cz um, und die Funken- 
strecke F3, die der Beleuchtung dient, schlagt 
durch. Die Verzigerung zwischen StobB- und 
Beleuchtungsfunken ist natiirlich von den Kapa- 
zitaten C, und Co, der Induktivitat L und den 


Abb. 3. Gesamtansicht des Gerates. 
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Elektrodenabstanden der Funkenstrecken ab- 
hangig. 

In Abb. 3 ist die Gesamtansicht des Gerates 
gezeigt. Vorn links die Hochspannungsquelle H, 
In der Bildmitte auf der optischen Schiene 
stehend die Wasserkiivette mit der Membran M, 
die aus mechanischen und optischen Griinden ein 
wenig gewolbt ist. Davor ein Projektionsobjektiv 
P. Beleuchtet wird die Kiivette durch die Fun- 
kenstrecke F3 mittels des Hohlspiegels S und des 
Kondensors K. Hinter der optischen Anordnung 
ist die Spule L zu sehen, die zum Zwecke der 
Verzégerungsregulierung mit mehreren Anzap- 
fungen versehen wurde. In das Innere der Spule 
lassen sich Transformatorbleche einschieben. 

In der Tabelle I sind die Daten der verwende- 
ten Teile und Gerate zusammengestellt. 


Tabelle I 


Spannung des Hochspannungsgenera- 
torsH . boxe) ic Sr eee ies 

KNapazitaten’ C775, Gat ee eee 

Abstand der Funkenelektroden F,, 


10 kV 
0,05, 0,5, 0,5 uF 


Eeations SA ALES. ta diet Pe ne eee 5, 0,2, 8 mm 
Dicke der Membran M (Uhrfeder). . | 0,4 mm 
Induktivitaét der Verzégerungsspule L | 2 mH---2 H 
Ladewiderstand R 50 kQ 


Kondensor K 


GW eae ems ello 
Projektionsobjektiv P . 


fi Ok cui li2 0 


4. Aufnahmen der Kavitationsblasen 


Im folgenden sollen einige Aufnahmen den 
Ablauf des Kavitationsvorganges, wie er oben 
geschildert wurde, illustrieren!. 

Fiir die erste Aufnahme, Abb. 4a, wurde die 
Verzégerungszeit zwischen StoB und Belichtung 
sehr kurz gewahlt, etwa 10-4s. Das Wasser hat 
sich eben von der Membran abgelést. An der 
Rauhigkeit seiner Oberflache erkennt man, dah 
das nicht ganz gleichmaBig geschehen ist, was 
auf die Wirksamkeit der eingangs erwahnten 
Keime an der Membranoberflache zuriickzu- 
fiihren ist. 

DaB nicht nur an der Membranoberfliche Kavi- 
tationskeime vorhanden sind, sondern auch in 
groBer Menge im Innern des Wassers, sieht man 
an den vielen schwarzen Punkten in der Bild- 
flache von Abb. 4a. Es sind kleine Kavitations- 
blaschen, die sich um die Keime gebildet haben, 
kurz bevor sich das Wasser von der Membran 
losriB; denn bis dahin war ja die ganze vor der 
Membran befindliche Wassermasse einer starken 


1 Natiirlich zeigt jedes Teilbild der Abb. 4 eine andere 
Blase. Es sind nicht kinematographische Aufnahmen einer 
einzigen. 
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Abb. 4. Aufnahmen von Kavitationsblasen bei verschiede- 
ner Verzégerung. 
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Zugbeanspruchung ausgesetzt. Nach dem Ab- 
reiSen der Oberflache 1a8t die Zugspannung so- 
fort nach, und die Blaschen implodieren augen- 
blicklich. Auf den spateren Aufnahmen treten sie 
in dieser Form und Gréfe nicht wieder auf. 

Die folgenden Bilder, Abb. 4b---g, sind mit 
stufenweise vergréBerter Beleuchtungsverzige- 
rung aufgenommen. In Abb. 4b ist die Ablésung 
des Wassers vollzogen, und in Abb. 4c beginnt 
bereits Wasser von den Seiten her nachzustrémen. 
Die weiteren Bilder zeigen den Zerfall der Blase. 
Ihre Oberflaiche ist merkwiirdig zerkliiftet, ein 
Umstand, auf den hier nicht eingegangen werden 
soll. Wahrend der letzten Zerfallsphase erscheint 
die Blase in einen Nebel bestehend aus vielen 
winzigen Blaschen gehiillt. Er ist auf den letzten 
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Bildern deutlich zu sehen und bleibt schlieBlich 
allein iibrig. Uber diesen Nebel, der eine typische 
Begleiterscheinung des Kavitationsvorganges ist, 
soll bei spaterer Gelegenheit noch berichtet wer- 
den. 


(Eingegangen am 20. Dezember 1954.) 
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ULTRASONIC ABSORPTION IN ANILINE-NITROBENZENE 
AND ANILINE-ETHYL ALCOHOL MIXTURES 


by D. SEeTrTE 


Istituto Nazionale di Ultracustica “O. M. Corbino”, Rome, Italy 


The ultrasonic absorption coefficient was measured in 
aniline-nitrobenzene and aniline-ethyl alcohol mixtures. In 
both cases a minimum of the coefficient was found at about 
0.3 aniline mole fraction. The results are shown in Fig. 1 
(x is the amplitude absorption coefficient and /f is the fre- 


quency). 


Alcohol 


—= (/F)10" 


Ethyl alcohol-aniline 
Nitrobenzene-aniline }—— 
t=25°C 


0 02 04 06 08 1,0 
—— Mole fraction aniline 
Fig. 1. 


In a previous work [1] a similar behaviour was found in 
mixtures formed by two unassociated highly absorbent li- 
quids having similar absorption coefficients, like benzene 
and carbon tetrachloride. It is possible in such a case to give 
a theoretical explanation of the experiments by assuming: 
1) that the dissipation in unassociated liquids is mainly 
due to relaxation in the distribution of energy between 
external and internal degrees of freedom of the molecules; 
2) that the probability of inelastic transitions in collisions 
of two molecules is larger when the molecules are of different 


kind, than when they are of the same type. It was shown [2], 
however, that this theory breaks down as soon as strong 
interactions appear either between molecules of any of the 
two components or between molecules of the two kinds. 

Aniline, nitrobenzene and ethyl alcohol are liquids in 
which a certain degree of molecular association is present 
and interactions among molecules of any two of them are 
to be expected. Therefore, the minimum found in the 
absorption coefficient curves cannot be explained in the 
same way as for mixtures of unassociated liquids. 

The temperature dependence of the absorption coeffi- 
cient was studied for a mixture aniline-nitrobenzene 
having 0.3 aniline mole fraction (Fig. 2). The definite 
decrease of the absorption coefficient shows that the dissi- 
pation is bound to some equilibrium destroyed by thermal 
agitation and confirms what has been said above. 

A similar dependence of the absorption coefficient on 
temperature has been observed in two categories of low 
viscosity liquids: 1) strongly associated liquids whose 
molecules can be found in various structural forms of 
different densities; it is the equilibrium among these 
structures which is affected by the pressure variations 
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produced by sound waves; 2) liquids in which a chemical 
equilibrium with a low rate constant is present; the equi- 
librium is this time altered by the temperature variations 
produced by sound waves. To the last category belong 
those binary mixtures [3] which have an absorption coeffi- 
cient curve with a maximum; the equilibria involved are 
those among different types of association of molecules of 
the two kinds. In case 1) the relaxation frequencies are 
well above the experimental range so that «/f? is practically 
constant with frequency; moreover Xexperimental ! X classical 
has a value roughly between 2 and 3 and it remains 
practically constant with temperature (for water [4] 
Xexp/%e is 3.11 at 20°C and 3.05 at 60°C; for methyl 
alcohol [5] the ratio is 2.9 at 20°C and 2.7 at 19.3°C). In 
case 2) the relaxation frequencies are usually much lower 
than in case 1) and o/f2 shows a definite decrease with f 
in the experimental frequency range. The ratio Xerp/X¢ 
has larger variations with temperature (for a water-methyl 
alcohol mixture with 0.25 alcohol mole fraction Xexp/x¢ is 
1.7 at 13.6°C and 2.7 at —10°C). 

On the basis of the temperature dependence of « in the 
aniline mixtures it seems plausible to admit that the excess 
absorption in these systems is due to relaxation phenomena 
of one of the two types mentioned. 

In order to have more indication on the kind of process 
taking place, we have calculated werp/x- at two tempera- 
tures for a mixture aniline-nitrobenzene with 0.3 aniline 
mole fraction and have measured «/f? as function of f. 

The values of density, viscosity and sound velocity 
needed for the calculation were obtained by means of 
direct measurements. The ration Xerp/x¢ is 2.5 at 18.5°C 
and 1.9 at —9.5°C. 

There is therefore a variation of Xerp/x¢- with tempera- 
ture; the average variation per degree (0.021), however, is 
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nearer to that found in pure methyl] alcohol (0.011) than 
to that observed in the water methyl alcohol mixtures 
above quoted (0.042). 

A study of absorption with frequency in the mixture 
aniline-nitrobenzene (0.3 aniline mole fraction) has given 
for (x/f2) x 101? the values 43.2; 43.3; 43.8 and 42 em7! -s® 
at the frequencies of 8; 25; 38; 53 Mc/s respectively. This 
shows that «/f2 is practically constant between 8 and 
53 Me/s. 

The two relationships Xerp/x¢ vs temperature and «/f? 
vs f seem to indicate that the relaxation phenomena which 
cause the excess sound absorption in the aniline-nitro- 
benzene mixture studied are mainly due to those structural 
equilibria which are altered by the pressure variations 
produced by sound waves. It seems therefore that an 
explanation for the behaviour of the absorption coefficient 
in aniline-nitrobenzene and aniline-ethyl alcohol mixtures 
must be sought on the basis of structural effects of this 
kind when either nitrobenzene or ethyl alcohol molecules 
are added to pure aniline. The strong hydrogen bonds 
which NHg groups of aniline can exchange with either NOz 
groups of nitrobenzene or OH groups of alcohol are pro- 
bably responsible for the origin of the structures formed 
by molecules of the two components. 

(Received 28th March, 1955.) 
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Bookreview 


L’acoustique des Orthoptéres — Compte rendu du 
ler Colloque International, Jouy-en-Josas (5—8 Avril 
1954), 444 pages, 17 x 24,5 em?, 1175 references, price 
fr. frs 2800. 


This book contains the papers read at an international 
symposium that took place in France in April 1954. Con- 
sequently the book gives a good review of current research 
in its special field. Being a collection of different papers 
by a large number of authors, some topics are discussed 
more than once, but this need not be a drawback. Unfor- 
tunately some of the contributions have to be considered 
only as preliminary reports so that not all are on the same 
standard. 

The papers are systematically arranged into three sec- 
tions. Section A deals with physical aspects of sound in gene- 
ral (units, definitions, ete.) and with the special topics which 
are of interest for the subject, especially the frequencies 
above 20 ke/s. Various types of apparatus for the recording 
analysis, and the production of these frequencies are de- 


scribed. The ‘tionophone” turned out to be a very valuable 
instrument for the production of high frequency sounds, 
but unfortunately it is nearly impossible to understand its 
construction and its action from the description given. 

Section B deals with anatomical details of the sound-pro- 
ducing structures. From the physical point of view it would 
have been of interest to have a similar careful survey of the 
sound-percepting organs. For a good understanding of the 
significance of communication by sound of insects, the 
study of the ‘‘soundfield”’ around a sound producing insect, 
also described in section B, is very important. 

Section C, finally, deals with the reactions of the animal 
on the sound produced by other insects of the same and 
other species and with the perception of these sounds. 

The book ends with a list of references of 44 pages! which 
should be almost a complete review of literature in this field. 
It will be evident from the foregoing review that the book 
is written for specialists; indeed they can find in it many 
valuable suggestions and information about new techniques. 


Hl. de Vries 


